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I. INTRODUCTION 


In this paper we will study the renormalization group (RG) flow of a simple Yukawa model 
describing relativistic fermions interacting through the exchange of scalar fluctuations. We will 
discuss some of its critical properties in a continuum of spacetime dimensions 2 < d < 4, 
dedicating most of the analysis to the d = 3 case. The class of models we want to consider is 
described by the following generic bare Lagrangian 

£ = ^ + U((/>) +^ V' ■ (I-l) 

where we have Nf copies of fermions, whose representation will be kept general in the following, 
and one real scalar field. The requirement of power-counting renormalizability would further 
restrict the interactions inside the potentials V and H (and would generically require the in¬ 
clusion of derivative interactions too) but we are not going to impose such conditions, since we 
are interested in describing the possible conformal models in this family, even if strongly cou¬ 
pled. In case the potentials V and H are even and odd respectively, the system is characterized 
by a chiral Z 2 symmetry, besides the \J{Nf) symmetry. For this reason, the model with bare 
potentials 

m = ^<(>'+ y, Hi<P) = y^ (1.2) 

is often called Gross-Neveu-Yukawa model, since it shares these symmetries with the purely 
fermionic Gross-Neveu model [T] and can be obtained from it by means of a Hubbard- 
Stratonovich transformation. 

Even for more general bare Lagrangians that are not related by any bosonization technique, 
the Yukawa models and chiral fermionic models remain deeply connected. The three dimensional 
Gross-Neveu model shows a second order quantum phase transition that separates the phase 
with preserved chiral symmetry from the one where this is spontaneously broken and a chiral 
condensate of fermions appears. The latter can be effectively described as a scalar degree 
of freedom, therefore this transition can be unveiled also as a dynamical effect in interacting 
scalar-spinor systems. Indeed, it is found that the critical properties of the Gross-Neveu model 
in 2 < d < 4 dimensions are compatible with the ones of the Yukawa model, thus indicating that 
the two are in the same universality class, which for generic but non-vanishing flavor number 
is also called the chiral Ising universality class. In both parameterizations, this is described 
by a non-Gaufiian fixed point (FP) of the RG flow. As a consequence, non-perturbative tools 
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are best suited for the investigation of its properties, and for the extraction of key quantities 
such as the corresponding critical exponents. Indeed several methods have been applied to this 
problem, including e-expansions lau, large Nj expansions lattice simulations [9HI3] 

and functional RG equations [HHiH]. 

These critical properties have great physical relevance for the description of several systems. 
In condensed matter, three dimensional relativistic fermionic systems, such as QED 3 and the 
Thirring model, play the role of building blocks for theories of high-Tc superconductivity [20) . 
and for the description of electrons in graphene |21j . Understanding the phase diagram and 
critical properties of these models at variable Nf represents pretty much the same challenge 
as the one posed by the Gross-Neveu and Yukawa model, and one can even address them in a 
unified picture [22] . Even the simple Yukawa model discussed in this work can hnd applications 
to extremely nontrivial phenomena in condensed matter. For the case of two massless Dirac 
fermions, its quantum critical phase transition in d = 3 might be a close relative of the putative 
transition between the semi-metallic and the Mott-insulating phases of electrons in graphene HE). 
For a single Dirac field instead, it is considered to be in the same universality class of spinless 
fermions on the honeycomb lattice with repulsive nearest neighbors interactions [13j . 

For a single Majorana spinor, it is a precious example of a three dimensional model showing 
emergent supersymmetry. Indeed, it is known that in this case the critical theory not only 
enjoys M = 1 supersymmetry, but also possesses only one relevant component, which means 
that by tuning a single macroscopic parameter one can discriminate between two distinct phases 
with preserved or spontaneously broken supersymmetry [I71|23|. On these grounds, a potential 
experimental realization of supersymmetry was proposed in |23| , at the boundary of topological 
superconductors. A similar phenomenon occurs for Yukawa systems with complex scalars and 
spinors, which have been argued to give rise to an emergent M = 2 supersymmetry |24j . 

The phase diagram of Gross-Neveu and Yukawa models has been analysed in d < 4 also for 
a better understanding of their d —?■ 4 limit. Glearly, nonperturbative phenomena in the latter 
case can have many applications in particle physics. These range from the chiral phase transition 
in QCD [25], where these models serve as simplified versions of quark-meson models [26], to the 
Higgs sector of the standard model [27] [28] , and to toy-models of composite-Higgs extensions [29] . 

In the present work we will analyze a more general truncation scheme for the functional 
renormalization group (FRG) study of these systems, showing under which conditions this brings 
important improvements in the results obtained by means of the latter nonperturbative method. 


3 


Such a truncation scheme amounts to allowing for a generic potential that essentially 

describes vertices with two fermions and an arbitrary number of scalars. This kind of interactions 
have been neglected in the FRG studies of fermionic models for a long time. Only recently they 
have been discussed in other works considering more complicated models and different but 
related questions. For example, in [30] the flow equations for this Yukawa system coupled to 
quantum gravity were derived, but only the linear coupling H ((/>) = ycj) was considered in explicit 
studies of these equations. Most prominently, in m the effect of higher Yukawa couplings on 
the chiral phase structure of QCD at finite temperature and chemical potential was analyzed by 
means of an effective quark-meson model. It was observed, within polynomial truncations of a 
Yukawa potential that higher order quark-meson interactions are quantitatively important 

in the description of the chiral transition. 

A similar but different study will be performed here, for the present Z 2 -symmetric Yukawa 
model, in lower dimensionality and for a generic number of flavors. We will confine ourselves to 
the study of the zero-temperature system at criticality, looking for scaling solutions for various 


d and Nf, and comparing the results obtained with different methods. In Sect. Ill we start with 
the leading order of the 1/iVj-expansion, reproducing known results in three dimensions, and 
generalizing them to multi-critical theories below three dimensions. Technical details regarding 
this analysis are sketched in App. [^ In Sect. IV we turn to a finite number of fermions and, 
by neglecting the wave function renormalization of the fields, we observe how critical Yukawa 
theories arise while continuously lowering the dimensionality towards two. To this end, we con¬ 
sider the FP equations for the two generic functions V{4>) and and solve them numerically 

without resorting to any truncation. In Sect. |Vj still neglecting the wave function renormal¬ 
izations, we adopt a different strategy for the numerical integration of the FP equations, and 
compute the global FP potentials in three dimensions, for various flavor numbers. For the case 
of a single Majorana spinor, we also apply these numerical methods to the computation of the 


critical exponents and perturbations. In Sect. VI we discuss polynomial truncations, showing 
how these can give results in satisfactory agreement with the global numerical analysis. As a 
consequence, we use them for a self-consistent inclusion of the wave function renormalizations, 
and produce estimates of the critical exponents in three dimensions and for various number 

we 


of fermions, which we compare with some of the existing literature. Finally, in Sect. VII 


address the d —>■ 4 limit at low number of fermions, and in Sect. VIII we draw a summary of 
our results. Yet, to introduce our work, we need to provide the reader with the definition of 
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the approximations involved in the computation of the flow equations, and with the resulting 
beta-functions. This is the object of the next Section and of App. 


II. THE RG FLOW OF A SIMPLE YUKAWA MODEL WITH 
MULTI-MESON-EXCHANGE. 


The functional renormalization group (FRG) is a representation of quantum dynamics based 
on Wilson’s idea of floating cutoff k. In this work we will adopt its formulation in terms of a 
scale-dependent IPI effective action, called average effective action [32]. For a given system, the 
form of this action is determined by the field content and by the symmetry properties, as well 
as by an initial condition (bare action) and boundary conditions for the integration of the flow 
equation 





(Ff)[$] + Rfc) 'kk 


(ii.i) 


Here (F^^^ [d>] represents the matrix of regularized propagators, while Rk is a momentum- 

dependent mass-like regulator. Since the dot stands for differentiation with respect to the RG 
time t = log k, this flow equation comprehends the inhnite set of beta functions for the infinitely 
many allowed interactions inside F^.. Extracting them amounts to projecting both sides of 
the equation on each separate interaction functional. In practical computations, one drops 
inhnitely many operators, thus performing a nonperturbative approximation called truncation 
of the theory space. To this end, several systematic strategies are available and appropriate in 
different circumstances, such as the vertex expansion or the derivate expansion. For reviews 
see [33] . 

In this work we will consider the following truncation: 


rk[(j),i>,-ip]= jd'^x + + . (II.2) 

Here cj) is a real scalar held, while V’ denotes Nj copies of a spinor held with real components. 
The latter parameter is related to the symmetries of the system and plays therefore a crucial 
role in pure fermionic as well as in fermion-boson models. Yet, as long as we truncate the 


theory space to the ansatz of Eq. (H.2), focusing on the mechanism of Z 2 -symmetry breaking, 
we can simply deal with the total number of real Grassmannian degrees of freedom Xf = d^Nf, 
considering it as an arbitrary real number. As soon as Xf >2 the truncation above is missing 
purely fermionic derivative-free interactions, that are indeed symmetry-sensitive and that would 
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contribute to the leading (zeroth) order of the derivative expansion. Furthermore, it is also 
missing field-dependent contributions to the wave functions renormalizations and which 
would appear in the next-to-leading (first) order of the derivative expansion. In the following we 


will call the ansatz of Eq. (II.2) a local potential approximation (LPA) for this simple Yukawa 
model, whenever the wave functions renormalizations are neglected = 1)) and 

therefore the fields have no anomalous dimensions = —dtlogZ^^^. The inclusion of the 
latter will be named LPAh Our justification for the choice of this truncation is in the exhaustive 
evidence that similar ansatze give a good description of the existence and properties of conformal 
models in 2 < d < 4 for linear systems with scalar degrees of freedom |33j . 


Projection of the Wetterich equation on the truncation of Eq. (II.2) yields the running of 


the corresponding parameters. Since we are interested in reproducing conformal models, that 
correspond to scaling solutions of the RG flow, it is useful to consider rescaled amplitudes 

kid-2)/2 }.{d-l)/2 


z 


1/2 




(jj 

since the new dimensionless renormalized field would then be constant at criticality. As a 
consequence we will focus on the potentials for these fields 

A/2 


Vk{4>) = k '^Vk 


Z" 


kid-2)/2 

In this new set of variables the flow equations read 


/c-i / zyv 


V = —c 


h = h{r]^-l) + 


d - 2+ r]^ 


= 




{(u(3))2 mf^\v") + 2Xf{h'f 


= 


d 



(II.3) 

)2/™(/l2,u")-/l"Zf^^(t'")} 

(11.4) 

i ) (pO 

(II.5) 

ph\v'')} 

^ 00 

(II.6) 


where Vd = {2'^^^7r'^^‘^T{d/2))~^, the threshold functions and on the right hand 

side denote regulator-dependent contributions from loops containing fermionic or bosonic prop¬ 
agators, and the equations for the anomalous dimensions are to be evaluated at the minimum (pQ 
of the scalar potential. Their definition can be found in App. together with the explicit form 
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they take for the linear regulator, which is our choice in this work since it allows for a simple 
analytic computation of such integrals. For this linear regulator the flow equations of the two 
potentials, read 


i, = -dv + -—+ Cd\ 
h = h{r}^-l) + -— ^ (/> h' 


1 - 


d+2 


1 + v'‘ 


-X 


1 - 


/■ 


■nj, ' 

d+l 


l + h 2 


(11.7) 


~\~Cd 


2h {h'Y 


1 _ 

^ d+1 


+ 


1 _ JX 

^ d+2 


(1 + ^2) (1+u") (1 +/l2) (1 + u") 


(i - Ife) 

(1 + v"Y 


(II.8) 


where we have denoted for convenience = Avd/d. 

A simple way of facilitating the stability of the vacuum is the requirement of Z 2 symmetry, i.e. 
invariance over </> —?• —cj). For standard Yukawa system, with a linear bare Yukawa interaction 
= ycj)^ this requires a discrete chiral symmetry tp ^ i'll) and tp -+ ip. A generalization of 
local interactions with such a symmetry then requires an odd There is also the possibility 

to let the spinors unchanged under the transformation, which would require an even function 
H{P). 

The goal of this work is to construct global FP solutions of the flow equations compatible 
with the symmetry conditions, and to study the properties of the RG flow in their neighborhood. 
The FPs, which describe scaling solutions, are computed by solving the coupled system of two 
ordinary differential equations u = 0 and /i = 0 or, in some cases, from the equivalent system 
for the quantities {v, y = h?). The dependence of such scaling solutions on the two parameters 
d and Xf is one of the main themes discussed in the literature as well as in the present work. 
Regarding the former, we will assume 2 < d < 4 and qualitatively discuss how the number of 
critical models varies with d, but we will especially concentrate on the properties of the d = 3 
system. For the latter, we restrict ourselves to non-negative number of degrees of freedom, and 
we start from the two simple limiting cases one can address. The simplest is Yj—)-0. In this 
case, the fermion sector remains nontrivial, see Eqs. ( II.4|II.6 ), but is not allowed to influence 
the scalar dynamics, which is therefore identical to the fermion-free model, see Eqs. ( II.3|II.5 ). 
Hence, as far as criticality is concerned, we expect to observe the same pattern of EPs that can 
be observed without fermions, with the same critical exponents in the scalar sector, even if at 
generically nonvanishing values of the Yukawa couplings. The second limit which brings radical 
simplihcations is Yj —)• 00 , and it is discussed in the next Section. 
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III. LEADING ORDER LARGE - Xf EXPANSION. 


Large-iVj methods are a traditional and successful way to analyze the strongly coupled do¬ 
main of the three dimensional Gross-Neveu model, which is renormalizable at any order in a 
l/A^/-expansion [HIT],[ 8 ]. As a consequence, any other nonperturbative method is challenged to 
reproduce known results in this limit. For this reason, before moving to the finite-Aj results 
provided by the FRG, let us start with discussing the behavior of this simple Yukawa model with 
many fermionic degrees of freedom, within the basic parameterization of its dynamics provided 


by Eq. (II.2), in a continuous set of dimensions 2 < d < 4. This FRG analysis, for the case of 
a linear Yukawa function, has already been performed in m- Our results can be considered as 
an extension of it, to include a generic function h{(p). As we will see, the main advantage that 
this brings at large-Aj is the possibility to describe also multi-critical models in d < 3. 

In this Section let us replace v with Xjv, as well as with Xfq^, and look at the leading 
order in 1/Aj. The first simplification is the fact that only canonical scaling terms and pure 
fermion loops survive. Therefore the flow equations at this order reduce to 


i! = -dv + - - ^ (pv' - 2vdl^P^{h?) (III.l) 

h = h{r,^-l) + h' (IIL2) 

(III.3) 

= 0. (III.4) 


Let us draw some general considerations about the FP solutions, by postponing the task of 
consistently solving the flow equation for The equation for h is almost regulator-independent 
(apart for the value of rj^j,) and the solution is a simple power 

h{4>) = Ch . (111.5) 


This is real only if the exponent is rational and with an odd denominator. Furthermore it is 
smooth only if the exponent is a positive integer. The FP solution for v is instead regulator 
dependent. Adopting the linear regulator, in 2 < d < 4 it reads 

v{4>) = c„ ,/,2'^/('^-2+’?0) _ ^ 2 F 1 (l, 1 - -h(0)2) . (IIL6) 

The function 2 ^ 1 ( 1 , 1 ~ —x), which actually can be reduced to a Hurwitz-Lerch function 

— (—X, 1, —0, has a logarithmic singularity at x = —1, therefore the condition that h{(f)) be 







real entails that this singularity is always avoided, and that the potential is globally defined. On 
the other hand, the smoothness of v is not for granted. Since 


2F1 



= 1 - 


d 


d 


d-2 


X — 


A-d 


X 


+ 0{x^ 


(III.7) 


and since this function is always convex, the leading ^-dependence of v at its minimum, i.e. at 
the origin, is provided by h'^{(j)) itself. Hence, the latter must be a smooth function, because 
we want the couplings associated to the derivatives of the potential at the minimum to be 
well dehned at the FP. The same reasoning, if applied to the Yukawa couplings, leads to the 
requirement that h((j)) be smooth at the origin. This translates into a quantization condition on 
the dimensionality of the scalar held 


d - 2 +rj^ 
2 


- , n G N 
n 


(III. 8 ) 


which is a consequence of the large-Xj limit. 

We hnd it helpful, for the interpretation of this relation, to consider a similar condition at 
the purely scalar FPs, with trivial Yukawa interaction. With this we mean the limit Xj ^ 00 
followed by Ch —>• 0 , which is clearly not the same as the fermion-free model; yet, by consistency, 
this limit should describe the classical properties of the latter model. Indeed, if Ch = 0 the only 
condition left is that the homogeneous part of the FP scalar potential be smooth and stable, 
that is 


d - 2 + _ d 

2 ’ “ ^ ’ 


n G N. 


(III.9) 


The meaning of this constraint is well known. By neglecting the quantum corrections, hence 
setting rj^ = 0 , one would deduce that the smooth bounded solutions u((/>) = a-re allowed 

only in 

2 r? 2 

dn = -- = 2 +-- , n G N . (III.IO) 

n — 1 n — 1 

This is the usual tree-level counting according to which the interaction cjp'^ is marginal in dn 
and becomes relevant for d < dn- From the quantum point of view, these dimensions are the 
corresponding upper critical dimensions for multi-critical universality classes. For any n, below 
dn a new FP with nontrivial r]^ branches from the Gaufiian FP and survives for 2 < d < 
dn |34t 135) . In the purely scalar model, this is already visible within a simple LPA of the FRG, 
where it is indeed possible to unveil and describe some properties of these universality classes in 
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a whole continuum of dimensions 2 < d < dm- In the leading order of the \aige-Xf expansion, 
the fact that quantum effects allow for these FPs at any 2 < d < dn remains invisible. This is 
because in the LPA one sets = 0, and in the LPA' the c/i —)• 0 limit again forces a vanishing 
anomalous dimension. This simply signals that the two limits Aj —>■ oo and h{(j)) —)• 0 do not 
commute. 

A similar analysis can be applied to the Yukawa system. Namely, if one forces classical scaling 
and sets 77 ,^ = 0 , the large-Aj limit constrains d to the critical values 

dn = 2 -\ — , uGN (III.11) 

n 

that are exactly the dimensions at which the interaction terms (/>" 7 /^ become marginal. Notice 
that they coincide with the critical dimensions of an even scalar potential, and that by selecting 
odd or even functions h{(j)) one can reduce the number of critical dimensions for h by a factor 
of two. As soon as anomalous scaling is allowed, the large-Aj limit tells us that the nontrivial 
FPs can indeed exist for d < dn, and quantizes the corresponding anomalous dimensions 


rjfj) = —\-2 — d = dn — d, n G N . 


n 


(III.12) 


Notice that they get smaller and smaller, the closer d is to the upper critical dimension dn- As a 
consequence, the value of Aj at which one expects a breakdown of the LPA with r]^ = 0 must be 
a decreasing function of (dn — d). Unfortunately, the latter is maximum for the very interesting 
n = 1 scaling solution, which includes the d = 3 Gross-Neveu universality class. However, even 
in this case, for small enough Xj we have no a-priori reason to discard the use of the LPA for a 
hrst study of the critical Yukawa models. On the other hand, for the n = 1 scaling solution the 
LPA' is able to reproduce Eq. ( |IIL12 ) and therefore provides a consistent picture of this critical 
model for any Xj, see App. This is not the case for the n > 1 multi-critical models, whose 
nontrivial scaling properties require larger truncations of the FRG. 

Before going on and discussing the finite-Aj results, let’s comment on the universal criti¬ 
cal exponents that one should approach in a large-A^ limit, since they provide an important 
reference point for the finite-Aj investigations. The eigenvalue problem for the linearized flow 
in vicinity of the large-Aj FPs is solved in App. both in the LPA and in the LPA'. The 
result is that one can safely split the problem into two classes of perturbations. The former have 
6h{<:f>) = 0 and Sv{4>) = 6cy(j)^, where we required the potential to be smooth, thus quantizing 
the corresponding critical exponents to the values 


9 m = d — M 


d-2 + r]^ 


= d- — 
n 


M gN 


(111.13) 
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i.e. the dimensionality of the couplings in front of 5v{(f)). The latter have 6h{(j)) = , a 

nontrivial 6v{4i), and again 


= jveN 


(III.14) 


where we used Eq. (III.8) as before. As a consequence, the large-Xj exponents are independent 
of Ch and c^. They are GauBian in the sense that they are directly linked to the dimensionality 
of the fields by naive dimensional counting, but the latter dimensionality, as far as the scalar is 
concerned, is deeply non-Gaufiian and actually independent of d. 

As usual one can observe a hierarchy among FPs with different n. For example, let us restrict 
ourselves to the slice of theory space parameterized by the couplings inside h{4>) only. Then, 
for a FP labelled by the integer n, there are n relevant operators, namely (fP,..., , and 

one marginal operator, itself. Within the LPA, the latter can be exactly marginal since it 
corresponds to shifts in Ch- For the n = 1 FP, the LPA' is enough to change this conclusion, since 
the flow equation for rj^ provides a condition that fixes the FP value of Ch- For n > 1, higher 
truncations are needed. Thus, the h-th FP can provide UV completion for theories approaching 
the n-th FP in the IR, only if n < h. The detailed study of the global flows among these FPs is 
in principle a straightforward task in the large-Aj approximation, but it is out of the purposes 
of the present work. We confine ourselves to sketching some properties of the FP potentials and 
of the linearized perturbations in vicinity of the FPs, which can be found in App. together 
with some comments on how these nontrivial critical theories disappear in d = 4. 


IV. LPA AT FINITE Xj AND GENERIC d. SOME FEATURES FROM NUMERICAL 


INVESTIGATIONS. 


In the previous Section we described how the large-Aj expansion supports the expectation 
that, as the number of dimensions is lowered from d = 4 towards d = 2, across the upper critical 


dimensions of Eq. (III.ll), new universality classes become accessible in the theory space of 


Yukawa models. In this Section we are going to present evidence that this happens also at finite 
Xf. Here and in the rest of this work, we restrict our analysis to the subset of theory space 
which enjoys a conventional Z 2 -symmetry, such that v is even and h is odd. Furthermore we 
adopt the LPA and neglect the flow equations for the wave function renormalization of the fields. 
As it was argued in the previous Section, as well as in App. with more details, one cannot 
expect this approximation to perform well for any n and Xj. Therefore the following studies 
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should be understood as a first step towards a proper description of these universality classes. 
Only the d = 3 chiral Ising universality class will be later analyzed also in the LPA', by resorting 
to polynomial truncations of the potentials, see Sect. m 

Since we look for odd Yukawa potentials, we can restrict the list of the operators that become 
relevant at the corresponding critical dimensions: 






dl{n > 2 ) 
dc{n > 0 ) 


2n _ 8 5 12 

“^’^’3’2’y ■■ 

4(n + 1) _ ^ 8 12 
2n + 1 3’y " 


(IV.l) 


In order to reveal the new universality classes appearing below these dimensions, we follow the 
strategy developed in [36], that was already successfully applied to the purely scalar model 
in continuous dimensions |35|. This consists in solving the FP condition, which is a Cauchy 
problem involving a system of two coupled second order ODEs, by a numerical shooting method, 
i.e. varying the initial conditions in a space of parameters which is two dimensional, since two of 
the four boundary conditions are fixed by the symmetry requirements (v'(0) = 0 and h(0) = 0). 
For the potential v we choose as parameter a = v”{0), relating it to u(0) using the differential 
equation. For h we use hi = /i'(0). Trying to numerically solve the non linear differential 
equations with generic initial conditions, one typically encounters a singularity at some value 
of 4>c{o',hi) where the algorithm stops. Such value increases in a steep way close to the initial 
conditions which correspond to a global solution, even if the numerical errors mask partially this 
behavior. As a consequence, in our case a three-dimensional plot for ^c(<7,/ii) is very useful to 
gain a first understanding of the positions of the possible FPs. 

In Fig. [^we show the results of this analysis, for Xj = 1 and for several dimensions: d = 
5,4, 3.9, 3.5, 3, I — For d = 5 and d = 4, as it is expected, we see a single spike in 

{a, hi) = (0,0) which corresponds to the Gaufiian solution. More details on this are given, for 


Xf < 1, in Sect VII In 3 < d < 4 we have crossed the threshold below which both the operators 
(j)'^ and become relevant, as is shown in Eq. (IV.l). In this interval, it is evident from the 
figure that we find three new spikes. One is characterized by di = 0 and cr < 0 and corresponds 
to the Ising critical solution. It is clearly visible in the fourth and fifth panels of Fig. but not 
in the third, since it is very close to the GauBian FP. The other two are physically equivalent, 
since they lie at opposite values of hi, and correspond to the chiral Ising universality class. They 
have (T < 0, which suggests that also these scaling solutions are in a broken regime for Xj = 1, 
at least in the LPA approximation. Moving to | < d < 3 we cross the marginality-threshold for 
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the operator but no other operators involving fermions have to be added to the set of the 
relevant ones. This corresponds to the appearance of the tricritical theory in the pure scalar 
sector, as we see from the new spike which develops with cr > 0 and hi = 0. Once d < | also 
the new operators cj)^ and become relevant and new critical solutions may appear. Indeed, 

in the left and the central plot of the third line of Fig. we see two new spikes, which again 
occur at opposite values of hi and are therefore equivalent, this time with cj > 0. Finally in 
the lower-right plot, where we present the case d = which is lower than | enough to clearly 
see the effects of the new relevant scalar operator 0®, one can appreciate the third new spike at 
u < 0 and hi = 0. The latter FP corresponds to the quadricritical scalar model as described 
for example in |35t 137). The former solutions, already assuming that they globally exist, define 
what one could call the chiral quadricritical Ising universality class, since they originate from 
the GauBian FP together with the purely scalar quadricritical model. 

We don’t show more plots with lower values of d, since the pattern is pretty clear. Pushing 
further this analysis towards dimensions close to d = 2, though conceptually straightforward, 
would probably anyway require more than the LPA. To provide the reader with some more 
details, in Fig. we zoom in the panel of Fig. that refers to d = | — . The three non trivial 

spikes which appeared at higher values of d > 3 are now out of this graph. From this figure 
one can see with more accuracy the presence of the three new nontrivial solutions. The two of 
them which lie at hi ^ 0, can also be visualized by a plot at constant value of u, approximately 
corresponding to the position of the peaks, see Fig. Here the range of hi is wider than in 
Fig.m so that one can see also a trace of the FPs generated at d < 4, which are nevertheless 
located at a different value of a. 

The analysis we discussed in this Section can be repeated for other values of Xj, thus getting 
a qualitative understanding of the position of the FPs as a function of both d and Xf. However, 
because of the uncertainties in the location of these peaks, it is hard to get a good qualitative 
knowledge of this function. Nevertheless, the latter is needed to prove that the arguments 
presented in this Section are rigorous, that each of the peaks corresponds to one FP, and to 
compute the corresponding critical exponents. For this reason, in the next Section we are going 
to adopt a different numerical method that will allow us to precisely answer these questions, 
focusing on d = 3 for definiteness, but allowing for a generic Xf. 
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FIG. 1: Spike plots for X/ = 1 on varying the dimension: d = 5,4, 3.9,3.5, 3, |, | — from left to 

right and from top to bottom. 


V. d = 3 LPA AT FINITE Xf. NUMERICAL SOLUTION OF THE FP EQUATIONS 


In this Section we construct, for some specific cases, the numerical solutions for v and h of the 


FP differential equations, obtained by setting Eqs. (II.7) and (II.8) equal to zero, in a domain 


for the dimensionless field (p that covers the asymptotic region. This is what might be called 
a global scaling solution. For convenience, we have actually considered the equivalent system 
for the quantities n(0) and y{(j)) = We focus here on d = 3 for which, from the analysis 

at Xf = 1 performed in the previous Section, we expect a FP with non-trivial scalar potential 
and Yukawa function. In the following we are going to take several values of Xf into account. 
After having found the corresponding nontrivial FP potentials, we determine the associated 
critical exponents and eigenperturbations. The knowledge of the global scaling solutions will 
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FIG. 2: Spike plot for d = | — and Xj = 1, zoomed area around the origin. 


<Pc 



FIG. 3: Spike plot for d = ^ ^ and Xf = 1, zoomed area around the origin. 


be important for a study of the quality of polynomial expansions, presented in Sect. VI . The 
latter approach is very useful especially in the case of the LPA', which gives us access to a 
self-consistent computation of the anomalous dimensions without enlarging the truncation to a 
full next-to-leading order of the derivative expansion. Clearly this programmatic analysis can 
be repeated for other values of d. 

We choose to construct a global numerical solution by starting from the knowledge of the 
asymptotic behavior allowed by the FP equations. Once the asymptotic expansions are deter¬ 
mined with sufficient accuracy we proceed, with a shooting method, to the numerical integration 
from the asymptotic region towards the origin. The properties of the solutions which reach the 
origin depend on the free parameters in the asymptotic expansions. By requiring the solutions to 
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transform correctly under Z 2 , one can uniquely fix the latter parameters to their FP values [35] 
The leading term of the asymptotic expansion for both v and h is determined, in the LPA with 
vanishing anomalous dimensions, by the classical scaling. Here we report the first correction to 
it. Denoting a = 2j{d — 2), the asymptotic behavior of the solution of the FP equations in the 
LPA reads 


?^asympt(</>) ^ A (j) 


2 q :-|-2 


2a ^d {B — 2AXf{a + l)(2a + 1)) 
^ 2AB{a + l){2a + l){d + 2) 

2a . i-2-2a Cd a{4:a{2a + 1)A + B) 


KsympM)- B(t) " + </> 


(V.l) 


+ 


2A2(a + l)(2a + l)2(d + 2) 
and depends on two real parameters A and B. In our analysis we have computed and used 
asymptotic expansions with eight terms for each potential. Starting the numerical evolution 
from some large value for (j) = ijimax, we have then investigated u'(0) and h{0) as functions 
of A and B. Computing numerically the gradient of these two functions, we were able to 
employ a kind of Newton-Raphson method to determine their zeros, i.e. the values of A and B 
corresponding to Z 2 -symmetric scaling solutions. In Fig. we present two examples of global 
solutions for the cases Xf = 1 and Xf = 2. The former is in the broken regime, since the 
Z 2 symmetric scalar potential has a non trivial minimum, while the latter is in the symmetric 
regime. Any solution (v, h) is characterized by two parameters, such as for example A and R, or 
v”{0) and h'{0), which indeed fix completely the Cauchy problem once they are complemented 
by the symmetry conditions. In Fig. we show the FP values of the integration constants A and 


B as defined by Eq. (V.l). The locus of the FP solutions in the plane (u'^(O), /i^(0)) as a function 
of Xf £ [10“^, 3] is instead presented in Fig. Notice that as Xf approaches zero, in the lower 
left end of the curve, h'{0) attains a hnite value, which is situated around 3.3. It is evident that 
the two regimes, broken and symmetric, are realized in two complementary intervals oi Xf. The 
transition between the two occurs at Xj ~ 1.64 for the LPA. In the next Section we will see that 
this value is slightly modified in the LPA', and becomes Xf cs 1.62. The vacuum expectation 
value (po and the value of h'{cfio) as functions oi Xf are presented in Fig. 

The critical exponents of these scaling solutions and the corresponding eigenperturbations 
are an important piece of information. This is obtained by studying the evolution of the small 
perturbations around the FPs. Therefore the linearized flow equations are the main tool to 
study such a problem. They are constructed, taking advantage of the separation of variables in 
4> and k, by substituting into the flow equations 

Vki(p) = + e6v{(j))e^^ , yk{(j)) = y*{(j)) + e6y{(j)) (V.2) 
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FIG. 4: The potentials v and h at the global scaling solution, computed numerically within the LPA. 
The case Xf = 1, which is in the broken regime, appears in the first two panels (top), while Xf = 2, in 
the symmetric regime, is shown in the last two panels (bottom). 



and then keeping the first term in e, for e <C 1. Such a procedure leads to the following eigenvalue 
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h'(0) 



FIG. 6: The values of (u"(0),/i'(0)) from the numerical global scaling solutions, varying in the range 
10“^ < Xf < One can notice the transition from the broken to the symmetric regime, which occurs 
at X/ ~ 1.64 for the present LPA. 



FIG. 7: The vacuum expectation value (j)o{Xf) from the numerical global scaling solutions is shown in 
the left panel, while in the right panel we plot the corresponding value of h'{(j)o){Xf), both in the LPA. 


problem 

0 = (A - d)6v + ^{d - 2)(i)5v' + Cd 
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and 


0 — (A — 2)5y + + Cd 


5v 


(2y {y + 1) - (y') ^ (y {v” + 5) + Sy^ + i)) 


y{l + yy{l + v"Y 


(3y2 + 2 y + 1 ) 


+'5y' y' 


(l + y)3K + l) 2 y 2 (i + y) 2 (i + ^//)2^ 
2 (3y + l) 


+ 


Sy'‘ 


{l+y)‘^{l+v") y (I y) (I + yllf ) {I + v"f 

where for simplicity we have renamed v* and y* as v and y. This system is of the form 


(V.4) 


(6 -a) 6f = 0. 


(V.5) 


if 6f is the vector of perturbations, 6f'^ = {6v,6y), and O is the corresponding differential 
operator. We have considered two different ways to solve this eigenvalue problem. 

The first approach is a direct generalization of the one we have already discussed for scaling 
solutions, in this case applied to the full set of equations: FP plus linearized flow. The asymptotic 
behavior of the eigenperturbations is computed by solving the asymptotic form of the linearized 
equations for large field, which is obtained using the known asymptotic expansion for v and y 


at the FP, given in Eq. (V.l). In d = 3 one finds 

<^^asympt 


6-2A ^ .-2A-4 + B-- (-2A^ + llA - 15)) ^ ^ 


135007rM2B2 


-8-2A 


(V.6) 


•^yasympt — 


^4_2a ^-2A-6 ( (2A2-11A + 15) (20^4 + B) ^ j3 (240.4A + B {2X‘^ + 5A-6))' 


168757r2yl3 


135007r2.425 


+0 


-10-2A 


In practice we used an asymptotic expansion with up to three terms per perturbation. We note 
that in a linear homogeneous problem the overall normalization of the eigenvector 5f plays no 
role. Therefore the asymptotic form of 6f depends only on a relative real parameter /?, which 
we choose to be a constant multiplying the leading term of 6y. One more free parameter is 
needed for tuning the behavior of the solutions at the origin, such that they fulfill the symmetry 
requirements (5u'(0) = 0 and Sy{0) = 0. This can be interpreted as the eigenvalue A itself. As a 
consequence, one expects a discrete spectrum of allowed values for A and j3. Unfortunately, due 
to numerical uncertainties, with this method we have been able only to restrict the eigenvalues 
to an interval described by a continuous function \{(5). Indeed one has to remember that the 
global numerical solutions have been constructed on some bounded neighborhood of the origin. 
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even if the latter overlaps with the region were the large field asymptotic behavior becomes 
dominant. Moreover, the linearized equations depend on derivatives of the numerical global FP 
solutions, for which the accuracy is reduced. 

The second approach we considered consists in inserting the known numerical FP solutions in 
the linearized equations, computing a numerical expression for all the (/)-dependent coefficients 
of this eigenvalue problem, and then solving them by means of a pseudo-spectral method based 
on Chebyshev polynomials. Also in this case some uncertainties remain, for the same reasons 
mentioned above. As an example, for Xf = 1 the leading critical exponent we find is 0i = — Ai = 
1.2279, which refers to the only relevant direction (we do not consider 6q = 3, since it is related to 
an additive constant in the potential and it is unphysical in flat space). All the other eigenvalues 
Aj are positive and associated to irrelevant operators, for instance O 2 = —A 2 = —0.6236 and 
03 = —A 3 = —1.5842. The relevant direction corresponds to the eigenperturbation 6fi = {6v, 6h) 
shown in Fig[^ Notice the fact that the relevant eigenpertubation has 6h[cj)) 7 ^ 0 unlike in the 
large-Aj analysis, where the only relevant perturbation compatible with symmetry requirements 
is 6v{(f>) = Scycj)'^, which corresponds to 9i = 1. Even if Xf = 1 is quite away from this 
limit, it is know that in this case the FP theory is a AA = 1 Wess-Zumino model [17( [23], and 
that the supersymmetry-preserving relevant perturbation is a change in the mass of the scalar 
field |17l I39j . which therefore leaves the Yukawa sector unchanged. Hence 5/i 7 ^ 0 is probably a 
consequence of the explicit breaking of supersymmetry introduced by our regularization scheme. 

We do not push further here the spectral analysis of the critical exponents and associated 
perturbations as a function of Xf, leaving it for a future study based on algorithms giving better 
control on the numerical errors. In the present work, these global numerical computations at 
Xf = 1 will serve as a reference for the development of a different, local, approximation method, 
based on polynomial truncations of the functions v{4>) and The latter will be discussed in 

the next Section, and will be also used for a more reliable discussion of the dependence of the 
critical exponents on the number of fermion degrees of freedom. 

VI. POLYNOMIAL ANALYSIS IN d = 3 

In this Section we are going to discuss the use of polynomial parameterizations and consequent 
truncations of the functions v{(f)) and h{4>). Though for definiteness we will address the specific 
case of the unique d = 3 nontrivial critical Yukawa theory, similar techniques can be applied 
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FIG. 8: Case d = 3 and Xf = 1: the components 6v and Sh of the relevant eigenperturbation, from the 
global numerical analysis of the LPA. 


to the other scaling solutions in 2 < d < 3, presumably with the same degree of success. 
Sect. VIA will present results obtained within the LPA, which can be directly compared to the 
full functional analysis developed in the previous Section. This will make us confident about the 
effectiveness and soundness of polynomial truncations, as well as of the necessity to go beyond a 
simple linear Yukawa coupling for an accurate description of critical properties of the theory. On 


these grounds, Sect. VIB will push forward the analysis to a self-consistent inclusion of the wave 
function renormalization of the fields, which is essential for quantitative estimates of the critical 
exponents, which will be compared with some literature for several values of Xj. Polynomial 
truncations will be also used in Sect. IVIII for some comments on the four-dimensional model. 

Let us start by presenting the truncation schemes we are going to analyze. Since we restrict 
ourselves to d = 3, we will demand v{4>) and h{4>) to be even and odd respectively. We will use 
the common notation p = ^^/ 2 , and we will adopt only one name for one and the same quantity, 
regarless of whether it is considered as a function of (/> or as a function of p. In the symmetric 
regime, the physically meaningful parameterization of the scalar potential is a Taylor expansion 
around vanishing field 

(VI.1) 


= E 

n=0 


Regarding the Yukawa potential, we are interested in two possible Taylor expansions, one for 
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h{(f)), already adopted in [3T], and one for y{p) = [h{(f))] . In the symmetric regime they read 


h{<p) = 


n=0 


Nh 


71=1 


yip) = 


(VI.2) 


(VI.3) 


In the regime of spontaneous symmetry breaking (SSB) the potential v{p) develops a nontrivial 
minimum k, = 0o/2, which becomes the preferred reference point for a different Taylor expansion 


iP) = ^0 + ^^{p-K) 


v[p] = 


(VI.4) 


n>2 


Though, in general, k is no special point for the function it still enters in the definition of 
the vertex functions, from which one extracts the physical multi-meson Yukawa couplings. As 
a consequence, in this regime it is necessary to change also the parameterizations of h{(f)) and 
y(/o), as follows 

Nh-l 


hr, 


h{4>) = 


n=0 


Nh 


71=1 


yip) = J2^_iip~ • 


(VI.5) 


(VI.6) 


The pair (Aly, N^), or more generally an ordering of the polynomial couplings by priority of 
inclusion in the truncations, can be chosen by relying on naive dimensional counting, as in an 
effective field theory setup, or on the knowledge of the dynamics at a deeper level, e.g. a global 
numerical solution for the FP functionals and the critical exponents. In the latter strategy one 
would sort the critical exponents in order of relevance and would try to accurately describe the 
corresponding perturbations. Alternatively, and maybe less efficiently, one could scan over the 
results produced by different pairs {Ny,Nh) and select them on the base of a comparison to 
the global numerical solution. In the former strategy instead, since the dimension of a scalar 
self-interaction is n, and the one of a multi-meson Yukawa coupling is 5/2 -|- n, 

we would expect that the pairs {Ny = D, = 17 — 2), for the truncation of h{(p) given in Eqs. 
( VI.2|VI correspond to including operators up to dimension D. However, since by truncating 
at level = D—2 we loose information about an operator of dimension Z7-I-1/2, if we want to be 
slightly more accurate we could include the latter and consider the pairs {Ny = D, = 17 — 1). 
In our analysis we did perform to some extent a random scan over different pairs {Ny, Nh), and 
we found that the two strategies nicely agree, so that {Ny = D,Nh = 17 — 1) is a very good 
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systematic choice for polynomial truncations. For similar reasons, as well as for the sake of 
comparison, we made the same choice also for the truncation of y{p) given in Eqs. ( VI.3|VlT6 ). 

It is necessary to stress that, in both the parameterizations given above, even at lowest order 
in the truncation for the Yukawa coupling, the beta-functions for ho or yi are different from the 
classic result [2^ illustrated in the reviews [33] and used for the present d = 3 critical theory 
for instance in [I1HI61[I81II9|. This happens because dth{4>), which comes from the projection 
of the r.h.s. of the flow equation on the term is a nonlinear function of (/>, independently 
of the parameterization of be it linear in (j) or not. Hence, in order to define the running 

of a linear Yukawa coupling, a further projection is needed. The prescription adopted by the 
above-mentioned studies is to identify the beta function of the linear Yukawa coupling with 
the first (/>-derivative of dth{^) at the minimum of the potential. For the truncations under 
consideration in this work instead, dtho comes from the zeroth order (/>-derivative of 
while dtyi is defined as the first order p-derivative of dty{p) = 2h{(p)dth{(^), always evaluated at 
the minimum of the potential. Simplicity is our main motivation for choosing a parameterization 
of the running Yukawa sector which does not include the traditional Yukawa beta-function, as 
we are now going to explain. 

The traditional projection has the structure of a Taylor expansion of dth{4>) about (f) = (po 
{(j)Q being the minimum of v{(p)). The choice of such an expansion for the parameterization of 
h{(f)) would entail an explicit breaking of Z 2 symmetry, which requires this function to be odd. 
Ideally, one would need to match two Taylor expansions, one about <p = (po and another one 
abont (j) = —(po, by imposing snitable conditions at the origin. These are just provided by Z 2 
symmetry. The result of this construction however is not a simple Taylor expansion any more 


h{(P) = '^ 0 )’" + (-1)"+^ {(p + (poY 

Zi Tl\ 

n=l 


(VI.7) 


and the projection rule on the generic coupling pn is more involved than simply taking the n-th 
0-derivative and evaluating it at 0 = 0o- Yet, it is true that the latter projection works for the 
Nfi-th coupling, such that this truncation does include the traditional beta-function of the linear 
Yukawa coupling as the Nf^ = 1 case. In this work we preferred to consider and compare only 
the two trnncation schemes presented in Eqs. (VI. 2|VI.5D and Eqs. ( |VI.3|VI leaving the one 


in Eq. (VI.7) aside. In the next Sections we are going to show that both polynomial truncations 
converge to the same results for large enough and N^, an observation that clearly should 
apply to all possible parameterizations. Fnrthermore, in both polynomial truncations simply 
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by setting = 1 one gets estimates that are significantly different from the full truncation- 


independent results. That the latter statement also applies to the truncation in Eq. (VI.7), can 
be assessed by comparison to the literature, which the reader can find in Sect. |VIB 


A. LPA 


In Sect.[V|we looked for the d = 3 nontrivial critical theories at varying Aj within the LPA, by 
means of numerical solvers for the ODEs defining the FP potentials. Here we repeat this analysis 
with the different method of polynomial truncations and we compare the results with the ones we 
previously found. The FPs emerge from the solution of a system of coupled nonlinear algebraic 
equations for the couplings. The critical exponents are defined by (minus) the eigenvalues of 
the stability matrix at the FP, i.e. the matrix of derivatives of the beta-functions with respect 
to the couplings [33]. The anomalous dimensions are computed in a non-self-consistent way, by 
neglecting them in the FP equations descending from Eqs. (II. and then by evaluating 

the flow equations for the wave function renormalizations Eqs. ( II.5[II.6 ) at this FP position. 

Let us start from the standard way of describing the Yukawa models, that is by approximating 
the Yukawa potential L((/)) with a single linear coupling. On the grounds of the results of the 
full functional analysis presented in Sect.jV] one could expect that this approximation performs 
well, since far enough from the large-held region the FP function h{(f)) does not strongly deviate 
from a straight line, see Fig. [^ For a linear Yukawa function, the expansions around the origin 
of h{4>) and y{p) give results which are identical order by order in Ny, both in the shape of 
the FP functions (in the sense that yi = 2/ig at the FP) and in the critical exponents. As a 
consequence we can present them in a single table for the former parameterization, the latter 
providing the same results. This is Tab. where we set e.g. Xf = 1. The hrst two critical 
exponents form a complex conjugate pair, which is clearly unsatisfactory. This is produced 
by the expansion around a trivial minimum of v{4>), that for Xf = 1 is not justified. Once 
we turn to the SSB parameterization of h{(j)), which is given on the left panel of Tab. |n| they 
become real. However, things become cumbersome for the single-coupling SSB parameterization 
of y{p), since we were not able to find any FP at all (which might nevertheless exist). Let us 
recall that, even in the case of a single Yukawa coupling, the beta functions descending from 
the two different polynomial truncations of /i((^) and y{p) are different, hence one cannot simply 
translate the FP position from one parameterization to the other. As soon as we add y 2 the FP 
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(N.,Nn) 

(2.1) 

(3.1) 

(-1.1) 

(5.1) 

(6.1) 

(7.1) 

(8.1) 

(9.1) 

(10.1) 


-0.04901 

-0.1225 

-0.1602 

-0.1743 

-0.1765 

-0.1740 

-0.1720 

-0.1716 

-0.1721 

h 

5.887 

6.841 

7.128 

7.204 

7.214 

7.203 

7.193 

7.191 

7.193 

h 

- 

84.22 

121.9 

134.7 

136.7 

134.5 

132.7 

132.4 

132.8 

h 

2.620 

2.464 

2.382 

2.351 

2.347 

2.352 

2.356 

2.357 

2.356 

Si 

1.701 

1.546 

1.438 

1.378 

1.358 

1.362 

1.372 

1.376 

1.375 

«2 

-1.050 

-1.156 

-1.246+i 0.2686 

-1.068+i 0.3386 

-0.9602+f 0.3238 

-0.9119+i 0.2933 

-0.9150+i 0.2844 

-0.9386+i 0.2941 

-0.9526+j 0.3044 

«3 

- 

-1.864 

-1.246-j 0.2686 

-1.068-i 0.3386 

-0.9602-i 0.3238 

-0.9119-j 0.2933 

-0.9150-i 0.2844 

-0.9386-i 0.2941 

-0.9526-i 0.3044 


0.2395 

0.2510 

0.2572 

0.2595 

0.2599 

0.2595 

0.2591 

0.2591 

0.2592 

J70 

0.2620 

0.2306 

0.2150 

0.2092 

0.2083 

0.2093 

0.2101 

0.2103 

0.2101 


TABLE I: Case d = 3 and Xf = polynomial expansion of around a trivial vacuum of the potential, 
with a fixed linear Yukawa function (standard Yukawa interaction), in the LPA. 


can be easily found. This then stimulates to consider the general effect of allowing for higher 
polynomial Yukawa couplings. 

The immediate observation is that their inclusion signihcantly alters the position of the 
FP and the critical exponents. Some degree of convergence is observed in several systematic 
strategies for the increase of and/or N^, but this can be convergence to the wrong results, 
i.e. to FP functions that do not agree with the numerical global solution. The linear Yukawa 
truncations provide one example of this fact. This is visible by comparing the two panels of 
Tab. [n| where on the r.h.s. we show the results provided by the {N^ = D,Nh = D — 1) 
systematic choice that we have already discussed above. The latter turns out to converge to 


the correct value of the FP couplings, as we are now going to argue. In Tab. Ill we show 
the results obtained by the systematic {D, D — l)-extension of polynomial truncations for y{p). 
Comparing the two panels one can see how the critical exponents can be computed by large 
polynomial truncations independently of whether these are around the origin or a nontrivial 
vacuum. Furthermore, comparing the right panels of Tab. Ill and Tab. |II] it can be observed 
how both the FP potentials and the critical exponents converge to values that are independent 
of the chosen parameterization. That these values are the ones corresponding to the full global 
solution provided in Sect. |Vj is shown in the right panel of Tab. Ill Notice however that there 
is a 0.6% difference between the relevant exponent computed with the polynomial truncations 
and the one obtained by the global numerical analysis. Even if we feel that we have the former 
method under a better control, we cannot give our preference to any of these estimates. 

In Fig. 1^ we plot different kinds of polynomial solutions, all in a (A/, = 9, = 8) truncation, 

against the numerical global FP functions, still for Xf = 1. For the potential v we show only 
the domain (p > 0.3, the agreement among all the curves being perfect for smaller values. The 
expansion around the origin has a smaller domain of validity as expected. Regarding the two 
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(iv„,jv,o 

(5,1) 

(6,1) 

(7,1) 

(8,1) 

(9,1) 

K 

0.01114 

0.01115 

0.01114 

0.01114 

0.01114 

^2 

25.08 

24.88 

24.80 

24.84 

24.85 

-^3 

813.8 

800.3 

793.33 

796.5 

797.5 

ho 

5.716 

5.690 

5.674 

5.681 

5.683 

Bi 

1.338 

1.333 

1.336 

1.336 

1.335 

62 

-0.2461 

-0.2466 

-0.2490 

-0.2484 

-0.2483 

93 

-2.232 

-2.060 

-2.033 

-2.067 

-2.075 

r]y, 

0.2629 

0.2288 

0.2288 

0.2288 

0.2288 

14‘ 

0.5259 

0.5166 

0.5155 

0.5160 

0.5162 


{N^,Nh) 

(5.4) 

(6,5) 

(7,6) 

(8.7) 

(9,8) 

K 

0.01002 

0.01009 

0.01008 

0.01007 

0.01007 

A 2 

15.34 

15.32 

15.30 

15.28 

15.28 

A 3 

508.3 

506.8 

503.6 

502.1 

502.1 

ho 

4.220 

4.211 

4.207 

4.206 

4.207 

h\ 

48.23 

47.73 

47.46 

47.43 

47.48 

0i 

1.231 

1.234 

1.236 

1.236 

1.235 

0-2 

-0.6144 

-0.6078 

-0.6080 

-0.6106 

-0.6117 

0.3 

-1.649 

-1.551 

-1.520 

-1.521 

-1.531 

'ny> 

0.3435 

0.3409 

0.3402 

0.3404 

0.3407 


0.4916 

0.4910 

0.4899 

0.4895 

0.4895 


TABLE II: Case d = 3 and Xj = 1, polynomial expansion of h{(j)) around a non trivial vacuum for both 
the potential and the Yukawa function, in the LPA, with or without the inclusion of multiple-meson- 
exchange interactions (right and left panel respectively). 



(5,4) 

(6,5) 

(7,6) 

(8,7) 

(9,8) 

( 00 , 00 ) 

K 

0.01000 

0.01013 

0.01006 

0.01006 

0.01007 

0.01007 

\2 

15.58 

15.17 

15.30 

15.28 

15.28 

15.28 

A 3 

521,8 

498.9 

503,0 

502.0 

502.3 

502.8 

yi 

44.59 

43.00 

43.51 

43.44 

43.43 

43,45 

yi 

1925 

1818 

1842 

1837 

1837 

1839 

01 

1.260 

1,221 

1.236 

1.236 

1.235 

1.228 

02 

-0.6849 

-0.7738 

-0.5964 

-0.6111 

-0.6127 

-0.624 

0.3 

-1.693 

-1.077 

-1.511 

-1.522 

-1.537 

-1.584 


0.3458 

0.3384 

0.3410 

0.3406 

0.3406 


V<t> 

0.4955 

0.4887 

0.4897 

0.4894 

0.4895 



(N^,Nh) 

(4.3) 

(5,4) 

(6,5) 

(8.7) 

(9,8) 

Ai 

-0.1209 

-0.1315 

-0.1339 

-0.1315 

-0.1309 

A 2 

10.60 

11.05 

11.16 

11.09 

11.06 

A 3 

293.2 

339.6 

351.0 

342.7 

340.1 

yi 

26.84 

28.38 

28.76 

28.53 

28.44 

y2 

986.6 

1161 

1206 

1178 

1167 

01 

1.324 

1.253 

1.226 

1.230 

1.236 

02 

-0.8293 

-0.7186 

-0.6410 

-0.5892 

-0.5989 

03 

-2.690 

-2.215 

-1.838 

-1.460 

— 1.446 


0.5209 

0.5615 

0.5716 

0.5642 

0.5618 


0.4486 

0.4645 

0.4683 

0.4663 

0.4654 


TABLE III: Case d = 3 and Xf = 1, polynomial expansion of y{p) in the LPA. Left panel: expansion 
around the origin, for which the global numerical solution provides Ai = —0.1313, yi = 28.47, and 
unstable higher couplings. Right panel: expansion around a nontrivial vacuum and, in the last column, 
the corresponding couplings extracted from the global numerical solution. 


set of expansions around a non trivial vacuum, the scalar potentials for the two cases are almost 
indistinguishable, while for the Yukawa function we obtain a slightly better result employing 


the one of Eq. (VI.6), as it is shown in the right panel of the hgure. The same kind of plots can 
be obtained for the polynomial truncations based on a single Yukawa coupling, corresponding 


to a linear Yukawa function. These are shown in Fig. 10, were we consider both polynomial 
expansions, around the origin and the non trivial minimum, for = 9. The left panel is 
especially interesting since it shows how, if one forces a linear Yukawa function, even with the 
SSB expansion, the shape of the potential is poorly reproduced. 

Having observed that in the LPA the {D, D — l)-systematic polynomial expansions converge 
to the global solution for Xf = 1, we assume that this is always the case, and make use of them 


for addressing how the FP and the critical exponents depend on Xf within the LPA. In Sect. Ill 


we have argued that when Yj is not small, there is no reason to trust the LPA for the d = 3 
critical theory, since rj^ should approch unity as Xf increases. This is what the global numerical 
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FIG. 9: Comparison of the Xf = I global numerical solution in the LPA (blue, continuous) with the 
corresponding {Ny = 9,Nh = 8) polynomial solutions, around the origin as in Eqs. ( VI.l| )-( VI.3) (red, 
dotted), around a non trivial vacuum as in Eqs. (VI.4 )-( |VI.6 ) (brown, dashed) and in Eqs. ( |VI.4 )-( VI.5) 
(green, dot-dashed), for the potential (left panel) and the Yukawa function y{(j)) = (right 

panel). 



EIG. 10: Comparison of the Xf = \ global numerical solution in the LPA (blue, continuous) with the 
corresponding {Ny =9,Nh = 1) polynomial solutions, around the origin as in Eqs. ( VI.l| )-( VI.3) (red, 
dotted) and around a non trivial vacuum as in Eqs. (VI.4)-( [VL5 ) (green, dot-dashed), for the potential 
v{(j)) (left panel) and the Yukawa function h{(j)) (right panel). 


analysis also indicates. Indeed in Sect. |V] we found that the constants A and B wildly grow 
from Xf = 3 on, in practice making the construnction of FP potentials harder and harder. This 
problem is easily addressed by means of the polynomial expansions. The results obtained with 
a (9,8)-truncation, both for h{(f>) and y{p), are shown in Tab. IV and Tab. [V| 

As expected, the anomalous dimensions show a very different Vj-dependence. Starting with 
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V, 

0.3 

0.6 

0.9 

1.2 

1.5 

1.64 

K 

2.311 10-2 

1.704 10-2 

1.173 10 -^ 

6.845 10-3 

2.219 10-3 

1.126 10-^ 

A2 

9.872 

12.21 

14.52 

16.75 

18.77 

19.61 

A3 

183.6 

294.3 

443.4 

632.5 

856.0 

967.6 

ho 

4.154 

4.178 

4,200 

4.218 

4.227 

4.230 

hi 

35,08 

40.29 

45.66 

51.12 

56.52 

59.04 

01 

1.435 

1.344 

1,261 

1.185 

1.117 

1.087 

O 2 

-0.6683 

-0.6481 

-0.6216 

-0.5896 

-0.5466 

-0.5212 

On 

-1.022 

-1.250 

-1.464 

-1.656 

-1.887 

-2.096 


0.2780 

0.3000 

0.3292 

0.3667 

0.4164 

0.4482 

V<t> 

0.2366 

0.3111 

0.4342 

0.6249 

0.8850 

1.014 


Xf 

1.64 

2 

2.5 

3 

3.5 

Ai 

-2.267 10-3 

0.1403 

0.5480 

1.705 

6.165 

A2 

19.50 

29.48 

65.58 

232.9 

1698 

A3 

960.5 

1955 

7265 

5.313 lO'^ 

1.090 10“ 

ho 

4.223 

4.600 

5.422 

7.041 

10.88 

hi 

58.84 

79.82 

142.6 

353.6 

1505 

0 i 

1.071 

0.9976 

0.9336 

0.9538 

1,041 

02 

-0.5212 

-0.4661 

-0.3727 

-0.2725 

-0.1783 

03 

-2.063 

-2.725 ±0.2953 

-2.763 ±0.8557 

-2.507 ± 1.242 

-1.956 ±1.695 


0.4521 

0.3372 

0.1066 

-0.1522 

-0.3048 

TI 4 , 

1.012 

1.545 

2.971 

6.660 

19.64 


TABLE IV: Case d = 3 and varying Xf, polynomial expansion of h((j}) around the non-trivial (left panel) 
or trivial (right panel) minimum for both the potential and the Yukawa function, with Nfi = 8 and Ny = 9 
in the LPA. 


X, 

0.3 

0.6 

0.9 

1.2 

1.5 

1.64 


Xf 1 1,64 

2 

2.5 

3 

3.5 

K 

2.310 in-2 

1.705 10"“ 

1.174 10-^ 

6.846 10-'^ 

2.187 10-^ 

3.115 10-® 


Ai 1-6.085 10-* 

0.1424 

0.5501 

1.706 

6.164 

A 2 

9.889 

12.21 

14.52 

16,75 

18.75 

19.56 


A 2 19.53 

29.52 

65.65 

232,8 

1698 

A 3 

184.1 

294.1 

443.4 

632.6 

853.3 

961.5 


A 3 1959.5 

1954 

7258 

5.301 10 ' 

1.089 10 “ 

yi 

48,27 

46.33 

44.26 

41,48 

37.82 

35.78 


yi 135.72 

42.37 

58.84 

99.13 

236.9 

y2 

1413 

1600 

1783 

1927 

1997 

1997 


y2 |l993 

2944 

6192 

1.990 10“ 

1.310 lO'"^ 

01 

1.436 

1.344 

1,261 

1.184 

1.112 

1.077 


81 |l.076 

1.003 

0.9374 

0.9551 

1.041 

02 

-0.6818 

-0.6643 

-0.6245 

-0.5897 

-0.5459 

-0.7877 


82 -0.5196 

-0.4652 

-0.3727 

-0.2726 

-0.1783 

03 

-1.021 

-1.242 

-1.467 

-1.665 

-1.864 

-0.5190 


03 '-2.006 

-2.582 

-2.794 ± 0.8023 

-2.520 ± 1.231 

-1.958 ±1.694 


0.2789 

0.2998 

0.3290 

0.3667 

0.4171 

0.4498 


7;i^|0.4509 

0.3360 

0.1061 

-0.1520 

-0.3048 


0.2367 

0.3111 

0.4342 

0.6249 

0.8850 

1.014 


Jl^|l.014 

1.548 

2.974 

6.659 

19.64 


TABLE V: Case d = 3 and varying Xf, polynomial expansion of y{p) around the non-trivial (left panel) 
or trivial (right panel) minimum for both the potential and the Yukawa function, with Nh = 8 and Ny = 9 
in the LPA. 


log yi 



FIG. 11: Behavior of the coupling yi in a Nh = 5, Ny = 6 polynomial truncation of y{p) around a trivial 
vacuum, within the LPA. The curve is a fit of data from Xf = 3.5 to Xf = 4 — 10“^. 


Vtp > for very small Xf, the former decreases and the latter increases as Xf is increased. 
Still for Xf around one, the two are small enough for qualitatively trusting the LPA, though for 
estimates of the critical exponents the LPA' provides different and more accurate results. The 
polynomial truncations agree with the global analysis and locate around Xf = 1.64 the transition 
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from the SSB to the SYM regime for the FP potential. Around this value rf^ reaches unity thus 
signalling the inconsistent use of the LPA. Yet if we insist on using this approximation for larger 
values of Xf, the breakdown of the approach is signalled by different phenomena. First of all 
the critical exponents become complex, from about Xf = 2 on. Then the anomalous dimensions 
rjff, and rj^, which are determined in a somehow un-legitimate way, become much bigger than 
unity and negative respectively. At the same time the couplings at the FP increase very rapidly, 
similarly to what was observed in Fig.[^ Actually in LPA it is easier than in the global numerical 
analysis to understand how quickly they grow. The result of a (6, 5)-polynomial truncation of 


y{p) around a trivial minimum is shown in Fig. 11 It is quite accurate to fit the behavior of 
the coupling yi close to Xf = 4 with a simple pole yi ~ 121.2/(3.999 — Xf). Also the remaining 
couplings have a rate of growth that is compatible to a divergence at a finite value of Xf, but 
these values would lie beyond the pole of yi. 

Also the comparison between the polynomial truncations and the global numerical results 
illustrates the appearance of severe problems as Xf increases. Moving to larger values ot Xf 
and entering the symmetric regime one sees, again comparing against the numerical solution of 
the ODEs, that the polynomial approximation has a smaller radius of convergence and therefore 
leads to a less trustworthy estimate of the LPA results. As an example we present the case 
Here the two curves show a good overlap for </> < 0.18, both for v{(p) and 


Xf = 2.5 in Fig. 


12 


y{4>), while at Xf = 1 the same grade of agreement was found for cj) < 0.28. Again the strongest 
restriction is imposed by the Yukawa function. Instead of interpreting these problems as a sign 
of the generic weakness of the polynomial truncations for large-Aj, we take the point of view 
that they are the way in which these truncations manifest the failure of the LPA for Xf roughly 
bigger than 1.6. We think that the results of the next Section support this interpretation. 


B. LPA' 

In the LPA' the anomalous dimensions are consistently determined by solving the FP equa¬ 
tions together with the flow equations for the wave function renormalizations. In the previous 
Sections we have shown that this is necessary for a correct qualitative description of the dy¬ 
namics of the model, roughly above Xf ~ 1.6. The expectation is that thanks to the wave 
functions renormalizations the system should gradually move towards the large-Aj limit, as it 
was already checked for truncations with a linear Yukawa function [HHIZ!. In this Section we 
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FIG. 12: Comparison of the numerical solution in the LPA (blue, continuous) with the corresponding 


(Ny = 9, Nh = 8)-polynomial solutions, for Xf = 2, around the origin as in Eqs. (VI.1)-(VL3) (red, 
dotted), around a non trivial vacuum as in Eqs. (VI.4)-(VI.6) (brown, dashed) and in Eqs. ( |VL4l )-( |VL^ 
(green, dot-dashed), for the potential v (left panel) and the Yukawa function y{4>) = (right panel). 



( 5 , 1 ) 

( 6 . 1 ) 

( 7 . 1 ) 

( 8 . 1 ) 

( 9 , 1 ) 

K 

6.250 10 -“ 

0.01261 

0.01262 

0.01262 

0.01262 

A2 

6.299 

6.995 

7.000 

7.001 

7.000 

A3 

52.38 

64.06 

64.28 

64.33 

64.29 

ho 

2.139 

2.533 

2.534 

2.534 

2.534 

Bi 

1.595 

1.548 

1.548 

1.548 

1.548 

S2 

- 0.7528 

- 0.6828 

- 0.6832 

- 0.6828 

- 0.6828 

B3 

- 1.241 

- 1.289 

- 1.299 

- 1.297 

- 1.294 


0.1168 

0.1273 

0.1272 

0.1272 

0.1272 

7]^ 

0.2807 

0.2237 

0.2238 

0.2238 

0.2238 



(5,4) 

(6,5) 

(7,6) 

(8,7) 

(9,8) 

K 

0.01080 

0.01078 

0.01077 

0.01078 

0.01078 

A 2 

6.009 

5.998 

5.997 

5.998 

5.999 

A 3 

61.01 

60.50 

60.47 

60.54 

60.56 

ho 

2.474 

2.473 

2.474 

2.474 

2.474 

hi 

7.548 

7.530 

7.542 

7.545 

7.544 

0 i 

1.444 

1.443 

1.443 

1.443 

1.443 

0-2 

-0.7721 

-0.7739 

-0.7745 

-0.7743 

-0.7741 

0.3 

-1.078 

-1.077 

-1.084 

-1.086 

-1.085 


0.1535 

0.1535 

0.1536 

0.1536 

0.1536 

7}d> 

0.2214 

0.2211 

0.2211 

0.2212 

0.2212 


TABLE VI: Case d = 3 and Xf = 1, polynomial expansion of around a non trivial vacuum for both 
the potential and the Yukawa function, in the LPA', with or without the inclusion of multiple-meson- 
exchange interactions (right and left panel respectively). 


want also to understand how big are the effects of the wave function renormalizations on the 
critical exponents, already for small Xf. 


As in the previous Section, let us start our discussion with the Xj = 1 model. Tab. VI is the 
LPA' version of Tab. |n| which considers the truncation of h{4>) with or without higher Yukawa 
couplings. If the effect of the inclusion of multi-meson exchange on the relevant exponent 0i was 
of the 8% in the LPA, it gets reduced to the 7% in the LPA'. However, in the truncation of y{p) 


the effect is of the 20%, see Tab. VII Also, the convergence of the polynomial truncations seems 


quicker in the LPA'. A comparison between the left panels of Tab. VI and Tab. VII illustrates 
how the predictions of the FRG can be made independent of the truncation scheme, here in 
the form of a different definition of Yukawa couplings, only by including full functions of field 
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{NyN„) 

(5,1) 

(6,1) 

(T,l) 

(8,1) 

(9,1) 

K 

9.208 10 -=' 

9.210 10 -" 

9.212 10 -" 

9.213 10 -" 

9.212 10 -" 

X2 

8.300 

8.307 

8.315 

8.316 

8.314 

As 

72.23 

72.45 

72.77 

72.82 

72.71 

yi 

18.64 

18.65 

18.67 

18.67 

18.67 


1.732 

1.731 

1.732 

1.732 

1.732 

«2 

- 0.5319 

- 0.5324 

- 0.5325 

- 0.5318 

- 0.5321 

«3 

- 1.626 

- 1.657 

- 1.676 

- 1.672 

- 1.664 

r]^ 

0.1886 

0.1887 

0.1887 

0.1887 

0.1887 


0.2680 

0.2681 

0.2683 

0.2684 

0.2683 


{N^,Nh) 

(5,4) 

(6,5) 

(7,6) 

(8.7) 

(9,8) 

K 

0.01079 

0.01077 

0.01078 

0.01078 

0.01078 

A2 

6.005 

5.997 

5.997 

5.999 

5.999 

A3 

60.83 

60.43 

60.50 

60.59 

60.56 

yi 

13.05 

13.04 

13.04 

13.04 

13.04 

yi 

152.0 

151.4 

151.7 

151.8 

151.7 

Oi 

1.444 

1.443 

1.443 

1.443 

1.443 

0-2 

-0.7710 

-0.7738 

-0.7745 

-0.7743 

-0.7741 


-1.072 

-1.077 

-1.086 

-1.086 

-1.084 


0.1536 

0.1536 

0.1536 

0.1536 

0.1536 


0.2214 

0.2211 

0.2211 

0.2212 

0.2212 


TABLE VII: Case d = 3 and Xf = 1, polynomial expansion of y{p) around a non trivial vacuum for both 
the potential and the Yukawa function, in the LPA', with or without the inclusion of multiple-meson- 
exchange interactions (right and left panel respectively). 


X, 

0.3 

0.6 

0.9 

1.2 

1.5 

1.62 

K 

2.377 in-2 

1 . 79310 -" 

1.253 10"^ 

7.316 10-3 

2.171 10-3 

1.164 10-^ 

A 2 

5.719 

6.028 

6.045 

5.849 

5.530 

5.385 

As 

55.00 

61.19 

61.55 

57.38 

50.81 

47.92 

ho 

2.745 

2.641 

2.518 

2.385 

2.252 

2.201 

hi 

9.355 

8.798 

7.890 

6.831 

5.789 

5.400 

01 

1.537 

1.490 

1.453 

1.427 

1.411 

1.407 

02 

-0.8158 

-0.7883 

-0.7755 

-0.7751 

-0.7833 

-0.7879 

0 ^ 

-0.9829 

-1.066 

-1.089 

-1.063 

-1.004 

-0.9742 


0.1510 

0.1529 

0.1537 

0.1531 

0.1514 

0.1505 


0.1366 

0.1687 

0.2073 

0.2499 

0.2936 

0.3108 


X/ 

1.62 

2 

3 

4 

6 

8 

Ai 

-7.622 10--^ 

4 . 13510 -" 

0.1443 

0.2316 

0.3602 

0.4448 

A 2 

5.375 

5.472 

5.604 

5.562 

5.185 

4.701 

As 

47.83 

43.65 

32.95 

23.64 

11.05 

4.560 

ho 

2.198 

2.157 

2.037 

1.915 

1.703 

1.538 

hi 

5.388 

4.863 

3.635 

2.694 

1.537 

0.9481 

01 

1.277 

1.229 

1.134 

1.077 

1.024 

1.004 

02 

-0.7776 

-0.7742 

-0.7794 

-0.7962 

-0.8345 

-0.8649 

03 

-0.8944 

-0.9581 

-1.101 

-1.196 

-1.287 

-1.311 


0.1508 

0.1314 

9.347 iir" 

6.939 icr" 

4.341 10-" 

3.073 10"" 


0.3106 

0.3721 

0.5057 

0.6024 

0.7223 

0.7894 


TABLE VIII: Case d = 3 and various Xf, polynomial expansion of h{cj)) around the non-trivial (left 
panel) or trivial (right panel) minimum for both the potential and the Yukawa function, with Nfi = 8 
and Ny = 9 in the LPA'. 


amplitudes, that is by allowing for higher polynomial couplings. 


Once we turn to the dependence of the results on Xf, which is shown in Tab. VIII and 


Tab. IX it becomes visible how the difference between the LPA and the LPA' can be negligible 
only for unphysical very small values of Xf. For 9i, it is the 7% at Xf = 0.3, and the 14% 


X/ 

1.62 

2 

3 

4 

6 

8 

Ai 

-7.366 10--* 

4.137 10-^ 

0.1443 

0.2316 

0.3602 

0.4448 

A 2 

5.374 

5.471 

5.604 

5.562 

5.185 

4.701 

As 

47.81 

43.63 

32.95 

23.64 

11.05 

4.560 

yi 

9.667 

9.304 

8.296 

7.338 

5.804 

4.733 

y2 

94.77 

83.91 

59.23 

41.28 

20.95 

11.67 


1.277 

1.229 

1.134 

1.077 

1.024 

1.004 

02 

-0.7775 

-0.7742 

-0.7794 

-0.7962 

-0.8345 

-0.8649 

03 

-0.8935 

-0.9578 

-1.101 

-1.196 

-1.287 

-1.311 

Vp 

0.1508 

0.1314 

9 . 34710 -" 

6 . 93910 -" 

4.341 10-" 

3.073 10-" 

VP 

0.3106 

0.3721 

0.5057 

0.6024 

0.7223 

0.7894 


X, 

0.3 

0.6 

0.9 

1.2 

1.5 

1.62 

K 

2 . 37710 -" 

1.793 10 -" 

1.253 10"^ 

7.315 10-’ 

2.169 10"^ 

1.125 10"^ 

A 2 

5.719 

6.028 

6.045 

5.849 

5.530 

5.384 

As 

55.00 

61.19 

61.55 

57.37 

50.79 

47.90 

yi 

17.51 

15.62 

13.67 

11.85 

10.26 

9.690 

y2 

214.7 

192.0 

162.1 

131.55 

104.5 

95.07 

01 

1.537 

1.490 

1.453 

1.427 

1.411 

1.407 

02 

-0.8152 

-0.7882 

-0.7755 

-0.7751 

-0.7831 

-0.7877 

03 

-0.9833 

-1.066 

-1.088 

-1.062 

-1.003 

-0.9727 

Vp 

0.1510 

0.1529 

0.1537 

0.1531 

0.1514 

0.1505 

VP 

0.1366 

0.1687 

0.2073 

0.2499 

0.2936 

0.3108 


TABLE IX: Case d = 3 and various Xf, polynomial expansion of y{p) around the non-trivial (left panel) 
or trivial (right panel) minimum for both the potential and the Yukawa function, with Nh = 8 and Ny = 9 
in the LPA'. 
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already at Xf = 1. On the contrary, as we will see later in this Section by comparing our results 
to the literature, the effect of the inclusion of higher Yukawa couplings decreases with incresing 
Xf. The transition between the SSB and the symmetric regime for the FP potential in the 
LPA' is around Xf = 1.62, while it occurs at Xf = 2.31 for truncations with a linear Yukawa 
function m- From these tables it also seems reasonable to expect that in the Xf ^ 0 limit 
the Yukawa couplings attain finite nonvanishing values, as it was observed already in the LPA, 
see Fig. Also, the trend in the change of 9i and rjf, is compatible with an approach to the 
corresponding Ising values, thus further supporting the discussion at the end of Sect. |IIj As far 
as the Yj- —7- oo limit is concerned instead, the smooth transition to the large-Yj exponents is 


evident in the right panels of Tab. VIII and Tab. IX 


Let’s now come to the comparison of our results with the literature. The classic methods for 
the investigation of the critical properties of the Gross-Neveu and Yukawa models are the e- and 
the l/Yj--expansions [SHE]. The former can be of great utility since both expansions around the 
upper and the lower critical dimensions give comparable results, such that d = 3 does not seem a 
too wild extrapolation. Yet, some treatment for these asymptotic series is needed. Resummation 
is unfortunately out of reach since they are known only up to the second or third order [HIS], 
apart for the anomalous dimensions for which the computations have been pushed up to the 
fourth order |6]. Polynomial interpolations of the two different e-expansions have been studied 
in [T8| for the case Xf = 8, and we report their results borrowing their notations, such that Pjj 
denotes a polynomial which is i-loop exact near the lower critical dimension, and j-loop exact 
near the upper. We also report the crude extrapolations that are obtained by simply setting 
e = I in the expansions of 6i = and Also the t/Nf expansion clearly needs some 

care, since we are interested in low number of fermions. Actually we are going to refer to this 
method only for Xf = 8 and Yj = 4, corresponding to Nf = 2 and Nf = 1 respectively. Again 
only the second or third order is known Eli. For the correlation-length exponent 9i = 


v ^ we 


adopt the Pade approximant used in |I8) . while for the anomalous dimensions we refer to the 
Pade - Borel treatment reported in [15] . 

The available FRG literature is rich and it offers a precious background on which we can 
measure the effects of the enlargement of the truncation discussed in this work. Essentially all 
the past studies considered the LPA', including a scalar potential and a simple linear Yukawa 


We made use of the formulas reported in [3], with typos corrected according to the observations of |18| . 
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coupling that can be considered as the first order in the truncation of Eq. ( |VI.7[ ) . The 

only exception in this sense is provided by the supersymmetry-preserving scheme that has been 
applied to the Xf = 1 case, which retained a full superpotential |39l HU HJ], thus including 
multi-meson exchange in the Yukawa sector, and sometimes was pushed to the next-to-next- 
to-leading order of the (supercovariant) derivative expansion. Also the choice of regulators is 
diverse, comprehending the linear, the sharp and the exponential ones (which in the tables 
we abbreviate with lin, sha, exp). In some studies the scalar potential was approximated by 
polynomial truncations in the symmetric regime, for which we provide the corresponding 
(Nw in case of truncations of the superpotential for supersymmetric flows). In others, that we 
label by Ny = oo (or Ny, = oo), the differential equations for the FP and the perturbations 
around it were solved by numerical methods, which are different from paper to paper. Our 
results are labeled by Nh > 1 . 

Other methods to which we can compare in special cases are Monte-Carlo simulations and 
the conformal bootstrap. Both of them can give high-precision computations of the critical 
exponents, but so far they have had a limited application to low-Xf Yukawa models. For 
Xf = 8 two lattice calculations of the critical exponents are available. One based on staggered 
fermions m, though ignoring a sign problem, provides results which are in good agreement 
with continuum methods, as it appears from Tab. An independent work applying the fermion 
bag approach [ 12 ], that is free from the sign problem, is instead offering very different results: 
u = 0.83(1), r /0 = 0.62(1), Tjfj = 0.38(1). In both cases it is not clear whether the symmetry of the 
lattice model is the expected one in the continuum limit Recently, another sign-problem-free 
simulation adopting the continuous time quantum Monte-Carlo method for a model of spinless 
fermions on a honeycomb lattice, provides estimates of the critical exponents of the chiral Ising 
universality class for Xf = A, i.e. a single Dirac field m- These results are compared to those 
emerging from the continuum methods in Tab. m Surprisingly they are much closer to our 


estimates for the case Xf = 2, see Tab. XII 


Regarding the latter case, notice that the results from m are affected by the absence of some 
terms in the flow equations that, being proportional to the vev of the scalar, become important 
for Xf < 2 Their effect significantly reduces the value of v. Since upon inclusion of multi¬ 
meson exchange the transition from the symmetric to the SSB regime occurs at lower values of 


^ We are grateful to H. Gies for informing us about these discussions. 
® See the discussion in m- 
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u 

fli 

Vcp 


FRG = 9, Nh = 8) lin (this work) 

1.004 

0.996 

0.789 

0.031 

FRG (]V„ = 3, Ni, = 1) exp [Tg 

1.016 

0.984 

0.786 

0.028 

FRG (]V„ = 6, iV,. = 1) sha [TU 

1.022 

0.978 

0.767 

0.033 

FRG (]V„ = ll.JVa = 1) lin [Tr] 

1.018 

0.982 

0.760 

0.032 

FRG (1V„ = oo.JVi. = 11 lin [TsI 

1.018 

0.982 

0.756 

0.032 

FRG (1V„ = oo.JVi. = 11 lin [T 9 I 

1.018 

0.982 

0.760 

0.032 

Monte-Carlo IIOII 

1.00(4) 

1.00(4) 

0.754(8) 

- 

l/Nf 2nd/3rd order [^1181 

1.040 

0.962 

0.776 

0.044 

(2 -1- e) 3rd order 

1.309 

0.764 

0.602 

0.081 

(4 — e) 2nd order 

0.948 

1.055 

0.695 

0.065 

P2.2 interpolated e-cxpansion 1181 

1.005 

0.995 

0.753 

0.034 

P 3.2 interpolated e-cxpansion 1181 

1.054 

0.949 

0.716 

0.041 


TABLE X: Critical exponents for Xf = 8. For a short description of the approximations involved in each 
method, see the main text. 



u 


V<!> 


FRG {N^ = 9, iV/, - 8) lin (this work) 

0.929 

1.077 

0.602 

0.069 

FRG {Nj, — 3,Nh — 1) exp I15II 

0.962 

1.040 

0.554 

0.067 

FRG (V„ = 00 . Va = 1) lin [T5][T9l 

0.927 

1.079 

0.525 

0.071 

Monte-Carlo 11311 

0.80(3) 

1.25(3) 

0.302(7) 

- 

l/Nf 2nd/3rd order [TllSlIiSl 

0.955 

1.361 

0.635 

0.105 

(4 — e) 2nd order 

0.862 

1.160 

0.502 

0.110 


TABLE XL Critical exponents for Xf = A. For a short description of the approximations involved in 
each method, see the main text. 


Xf, our computations are still in the symmetric regime. This might qualitatively explain the 
drastic departure from the results of m- 


Also the comparison for Xf = 1, which is presented in Tab. XIII requires some comments. 
Let us recall that for this held-content the system at criticality is described by a A/" = 1 Wess- 
Zumino model [I71E3]. Hence, if the regularization does not break supersymmetry, the critical 
anomalous dimensions of the scalar and of the spinor should be equal. Furthermore, a superscal¬ 
ing relation = {d—T])/2, which was first observed in [JO] and later proved to hold at any order 
in the supercovariant derivative expansion in m. is expected to hold. This is what happens for 



u 


V4> 


FRG = 9,Nh = 8) lin (this work) 

0.814 

1.229 

0.372 

0.131 

FRG {JV„ = 3,Nh = 1) exp [Tsl 

0.633 

1.580 

0.319 

0.113 

FRG (JV„ = 3, JVfc = 1) lin [Ts] 

0.623 

1.605 

0.308 

0.112 

FRG {Nj, = 00 , Nh = 1) exp 1151 

0.640 

1.563 

0.319 

0.114 

FRG (JV„ = 00 , JVa = 1) lin [H] 

0.621 

1.610 

0.308 

0.112 

FRG (Ar„ = 00 , JVa = 1) lin [To] 

0.4836 

2.068 

0.3227 

0.1204 

(4 — e) 2nd order 

0.773 

1.293 

0.317 

0.154 


TABLE XII: Critical exponents for Xj = 2. For a short description of the approximations involved in 
each method, see the main text. 
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v 

Si 

S2 

770 

770 

3-201 

FRG (7V„ = 9,Nh = 8) lin (this work) 

0.693 

1.443 

-0.796 

0.154 

0.221 

0.114 

SUSY FRG (7V„ = oo) opt ji = 2 NLO gT] 

0.711 

1.407 

-0.771 

0.186 

0.186 

0.186 

SUSY FRG (7V„ = oo) opt ji = 2 NNLO gT] 

0.710 

1.410 

-0.715 

0.180 

0.180 

0.180 

SUSY FRG (7V„ = oo) opt ji = 1 g2] 

0.708 

1.413 

-0.381 

0.174 

0.174 

0.174 

SUSY FRG (7V„ = oo) opt n = 2 g2] 

0.706 

1.417 

-0.377 

0.167 

0.167 

0.167 

FRG (N„ = 2, TV,. = 1) 1-loop [Tt] 

0.72 

1.39 

-0.71 

0.15 

0.15 

0.22 

(4 — e) 1st order 1231 

— 

— 

— 

0.143 

0.143 

— 

(4 — e) 2nd order 

0.710 

1.408 

— 

0.184 

0.184 

0.184 

Conformal Bootstrap 1431 

— 

— 

— 

0.13 

0.13 

— 


TABLE XIII: Critical exponents for X/ = 1. About the FRG results, the schemes, the regulators, and 
the approximations are very different, see the main text. 


example in the e-expansions or in the SUSY FRG. Since the scheme adopted in the present work 
explicitely breaks supersymmetry, we expect and we observe violations of these properties. Also 
in HZ! supersymmetry is broken by regularization, and these violations are present, but they 
could be partially reduced or canceled by tuning the regulator. This tuning gives the results 
reported in Tab. pain A similar analysis of the regulator dependence of universal quantities 
and of the consequent breaking of supersymmetry could be performed in future studies for the 
present family of truncations. Yet, even by explicitly breaking the FP supersymmetry, we get 
exponents which are not very far from the ones produced by the above mentioned methods. Let 
us add few details on the SUSY FRG results shown in Tab. Xlllj They are obtained by setting 
one of the regulators to zero, and choosing a shape similar to the linear regulator for the other, 
with an exponent n that differentiates between the conventional linear regulator (opt n = 2) and 
a slight variant (opt n = 1). Also the truncation scheme is different from the one discussed in 
the present paper, since it is related to an expansion in powers of the supercovariant derivative, 
that has been considered at the level of the LPA' |39l I42| . at next-to-leading order (NLO) or at 
next-to-next-to-leading order (NNLO) [41j . For the case Xf = 1 we can also compare with a 
pioneering study based on the conformal bootstrap |l3|. In Tab. jXIII we included the one-loop 
computations of [HIES], even if two-loop results are on the market |3| , on the base of the naive 
observation that for Yukawa systems with complex scalars and spinors, whose FP should effec¬ 
tively show M = 2 supersymmetry [ 23 ] , the anomalous dimensions obtained from the first-order 
of the (4 — e) expansion, = 1/3, agree with the available exact results 
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FIG. 13: Spike plots for Xj = 0, v{(j)) = 0 and d S {3.5,3.7,3.9,3.99,3.999} from red (upper) to blue 
(lower) in the LPA'. 


VII. D=4 


From the leading order of the 1/X ^-expansion one expects that for large enough Xf the chiral 
Ising FP merges with the GauBian FP in the d —4 limit. Also at Xf = 0, for which we know 
from the discussion at the end of Sect. that only mirrored images of the purely scalar FPs can 
exist, one can observe that the latter merge with the GauBian FP for d —?• 4, compatibly with 


the presumed triviality of four-dimensional scalar theory. This is illustrated in Fig. 13, which 
is produced as Fig. but integrating only the FP equation for h{cj)) at v{cj)) = 0 and Xf = 0 
in the LPA'. Yet, it remains to be shown what happens for a small non-vanishing number of 
fermions. Dimensional analysis indicates d = 4 as the upper critical dimension for any Xf. This 
can be checked by means of the FRG, either by numerical integration of the FP equation, as 


it was shown for example in Sect. IV for Xf = 1, or by the polynomial truncations discussed 
in the last Sections. Indeed, the latter have already been used in the past, precisely to address 
this question. 

In fact, an exploratory study of what happens to the d —?• 4 limit in a Z 2 -symmetric Yukawa 
model with very small Xf was performed in m, in order to test a mechanism for the generation 
of nontrivial FPs in fermion-boson models, that has subsequently found in chiral-Yukawa models 
some natural candidates [28]. That analysis pointed out that within a (Y^ = 2,Y/j = I) 


polynomial truncation, according to the scheme of Eq. (VI.7), the FRG detects nontrivial FPs 


also in d = 4, for unphysical small values ot Xf. This holds both in the LPA and in the LPA'. 
However, the fact that the FP position and the critical exponents are significantly different in 
the two approximations was interpreted as a signal of the need to include further boson-fermion 
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interactions in the truncation, in order to understand if these FPs are physical or merely an 
artifact of the approximations. This Section reports on the changes brought by the different 
treatment of the Yukawa sector presented in this work. 

At the level of the LPA we generated three-dimensional plots similar to the ones illustrated 
in Fig. second panel, by shooting from the origin with random values of {v"(0), h'(0)), for 
several values of Yj < 1, and we looked for spikes signaling possible FPs, but we have not found 
any of them. We were also not able to produce any global solution studying numerically the 
Cauchy problem from the asymptotic region, along the lines of Sect. |Vj We then re-considered 
the analysis at the level of polynomial truncations. Already trying to reproduce the results 
of [27j in other truncations with Ny = 2 and Nh = 1, can be a nontrivial test, because of the 
different beta-function of the Yukawa coupling, associated to different projection rules. We have 
already argued that a change of the results depending on the parameterization employed signals 
the presence of errors induced by the use of inconsistent truncations. We first concentrated on 
the LPA, which at least for d < 4 is able to reproduce the right number of nontrivial FPs. In 
this case, the truncation adopted in m allows for non-Gau6ian FPs approximately for Xj < 1. 
For instance, at Xf = 0.4 one can find the following FP 

K = 0.00165 , A 2 = 27.26 , gf = 81.13 (VII.l) 

with two relevant directions 


01 = 2.372 , 02 = 0.592 , 03 = -2.859 . 


(VII.2) 


We observed that in a polynomial truncation of y{p) as in Eq. (VI.6), the FP position is different 


K = 0.00167 , A 2 = 54.18 , yi = 494.0 


(VII.3) 


as well as the critical exponents 


01 = 1.653 , 02 = 0.932 , 03 = -3.445 . 


(VII.4) 


Still, the changes are not dramatic. On the other hand, we could not find any real FP within the 


same order of the truncation of h{4>) given in Eq. (VI.5). We tried to circumvent this problem 


as in d = 3, by following the FP found in one parameterization to higher orders, and then 
translating back to the other parameterization. Yet, we were not able to reveal the FP for y{p) 
for bigger values of and Nh, nor to find it by chance in different orders of the truncation of 
h{(p). 
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Hoping that the inclusion of the wave function renormalizations could stabilize the polynomial 
truncations and help us in the search for FPs, we then considered LPA', using the results 
of [27] as a guide for the localization of the interesting region in the space of couplings. While 


the FP is present in the first order of the truncation of Eq. (VL7), we could not find it in 
the parameterizations considered in this paper. Let us once more stress that this does not 
completely exclude that it can be found by other methods, even if we consider this very unlikely. 
Nevertheless, for the LPA' we have not tried a numerical shooting at nonvanishing Xf as in 
the LPA. Hence, a more careful numerical analysis is needed, to exclude with a higher level 
of confidence the presence of low-Xf FPs in the theory space described by the truncation in 


Eq. (H.2). A even better test would be to consider the full next-to-leading order of the derivative 
expansion. 


VIII. CONCLUSIONS 

A proper quantitative control of the quantum dynamics of the Z 2 -symmetric Yukawa model, 
beyond the domain of applicability of perturbative methods, is important not only from a 
generic field-theoretical point of view, but also for phenomenological reasons, since the latter is 
very useful as a toy-model of numerous condensed matter systems, as well as of specific sectors of 
modern particle theory, see Sect.|l]for more details. The functional renormalization group (FRG) 
is a simple analytic nonperturbative method that can provide a detailed description of strongly 
coupled systems, under approximations that are testable and improvable in several systematic 
ways. Furthermore, these results can be produced, almost simultaneously, in a continuous 
number of spacetime dimensions d and fermionic degrees of freedom Xf, thus allowing for a 
quick analysis of the dependence of the dynamics on the latter parameters. 

In this work we focused on the critical behavior of the Z 2 -symmetric Yukawa model at zero 
temperature and density. Our principal aim was to test the impact of multi-meson-exchange, 
encoded in a Yukawa coupling which is a full function of the scalar field, on the FRG description 
of the latter behavior, a question that to our knowledge has never been considered before. 
Nevertheless, our analysis is relevant not only for the FRG community. For instance, in Sect. m 
we discussed the leading order of the 1/Yj-expansion, whose results can be directly exported out 
of the FRG framework in which we produced them, and recovered also by different methods. This 
study illustrated how, by allowing for multi-meson-exchange, one can describe the generation of 


38 




multi-critical conformal Yukawa models as d is lowered from d = 3 towards d = 2, across the 
corresponding upper critical dimensions dn = 2 + 2jn, with n a positive integer. We also showed 
how the large-Yj limit quantizes the corresponding critical anomalous dimension = dn — d. 


In Sect. IV we checked that this pattern of generation of critical theories as a function of d 
holds also at Xf = 1, and presumably at any other finite Xf. This would imply that in the 
d —7- 4 limit only the Gaufiian fixed point (FP) suvives. The latter statement, being of special 


relevance for particle physics, was further analyzed in Sect. VIT where we argued that it applies 


also for Xf < 1, at least within the ansatz of Eq. (II.2). Let us remark that, as far as we know, 
the observation of multi-critical conformal Yukawa models at finite Xf and in continuous fractal 
dimensions 2 < d < 3 is a novel result. 

Concerning the finite-Yj results, they indicate that in several cases the effect of multi-meson- 
exchange cannot be neglected, either quantitatively or even qualitatively. We argued that these 
higher Yukawa interactions are required by consistency of the truncation, otherwise the solutions 
of the system of differential equations defining the flow of the scalar potential v{(j)) and of the 
Yukawa “potential” would depend on the chosen parametrization of these functions. For 
instance, the same FP solutions should be reproduced using any polynomial truncation of these 
functions, at least within a certain domain. On these grounds we believe that general FRG 
studies of Yukawa models should at least consider the inclusion of these interactions, and possibly 
check when they can actually be neglected. 

On the quantitative side, in Sect. |V]we explicitly numerically constructed these global FP 
solutions for d = 3 and several values of Xj. These results include the Gross-Neveu univer¬ 
sality class for Xf > 1, and the superconformal J\f = 1 Wess-Zumino model for Xf = 1. At 
Xf = 1, we also numerically computed the critical exponent v, and the corresponding linear 


perturbation around the FP. In Sect. VI we showed how the results of the global analysis can be 
easily reproduced by two different kinds of high-order polynomial truncations. However, these 
studies were performed in the local-potential approximation (LPA), that is by neglecting the 
renormalization of the fields. Taking into account the anomalous dimensions (LPA') was crucial 


to obtain a more accurate picture, especially for Yj > 1, so that in Sect. VIB we developed a 
LPA' analysis, based on the same polynomial truncations which were proved to be trustworthy 
in the LPA studies. 

This allowed us to to produce estimates of the critical exponents and in d = 3 and 
various Xf, and to compare them with some of the existing literature. We concentrated on the 


39 







especially interesting cases of two and one massless Dirac {Xf = 8 and 4), of a Weyl {Xj = 2), 
and of a Majorana spinor {Xf = 1). They can be found in Tab. |X, XI, XII XIII Often, there 
still exists some significant mismatch among the available estimates, such that more studies by 
all kinds of methods, including Monte-Carlo simulations and higher-order e- or 1/Xj-expansions, 
are welcome. As far as the FRG is concerned, the results seem stable for Xf > A, while for 
lower number of fermions there is still room for debate, and probably larger truncations are 
needed. The supersymmetric case Aj- = 1 is an exception also in this sense, since it enjoys a 
good agreement among the results produced with different methods. 

Larger truncations, such as a next-to-leading order of the derivative expansion, are anyway 
needed for a quantitative analysis of multi-critical models in 2 < d < 3, as we argued in 


App. Ill in the large-Aj limit. Still within the LPA', the next natural step is to produce global 
numerical studies similar to the ones presented for the LPA in Sect. IV and |Vj Regarding the 
possible applications of the present analysis to different models, one possibility is to enlarge the 
symmetry group from Z 2 to 0(N). The N= 3 three-dimensional chiral Heisenberg universality 
class, for instance, can be interesting for the physics of electrons in graphene [18]. With an 
enlarged symmetry, the effect of different representations would become a natural case-study 
and would further widen the class of physical applications of these studies |1S| . The same kind 
of truncation can also be used in the context of a Yukawa model interacting with gravity, along 
the lines of [l6], to investigate first the asymptotic safety properties of the model, and then to 
construct global flows from the UV to the IR. Some scenarios could be of particular interest for 
cosmology. 
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Appendix A: Regulators and threshold functions 


A -K (2) 

We have to evaluate the r.h.s of Eq. (II.1), for which we need the matrix. Considering 

the held ip as a column and ip as a row vector, let us denote by the row vector with 

components (P{q)i arid by d>(p) the column vector given by its transposition. 

( 2 ) 

Then T^ ^ is obtained by the formula 


p(2) _ 

k — 


:^k- 


This inverse propagator is regularized by addition of the following regulator 

Rb{p) 


Rk{q,p) = S{p-q) 


0 



where 


Rb{p) = Z^p’^rjiip^), 

( 0 (5dyT\ 

-Rf(p) = - Z^rpip^) , 

0 / 

is a 2d^Nf x 2djNf matrix. In principle one can have different regulators for the scalar (B) and 
for the spinors (F). A compact way to rewrite the how equation is 

9irfc = ^4sT:iog(rf+ d?fc), 

where 

~ ^ dtiZ^rp) 6 dtjZ^rp) 5 

‘ ” Z^ drp Z^ 5rp 

and • denotes multiplication as well as integration over the common argument of the shape 
functions of the two factors. Then the regularized kinetic (or squared kinetic) terms are given 
by Pbip{x) = x(l + rB/F(x)), and the loop momentum integrals appearing on the r.h.s. of the 
how equation give rise to corresponding regulator dependent threshold functions. Introducing 


41 







the abbreviation fp^f these threshold functions read 

^ J dt log (Pb/f + 




„2-d 


4fd 


dt 


Ip -Pb/f 


+ ojk?' 




‘^Vd 


(Pp + oJik^Y^ (Pb + 0J2k‘^Y'^ 


II 

k6-d 

to 

‘^Vd , 

II 

ki-d 

[ P^dt 

Jp 


{ooi) = 

^6-d 

to 

4rd , 


k^-d 

to 

4ud , 


d 


dp^ Pp + Lvk"^ ^ 
4a / ^ 1+rp y 


9p' 


rPB 


(Pb+wiP)^^ 

1 + rp 


wP^ 


Pp+Wife2 (Pb+W2A;2)2 


In this work we adopted the linear regulator xrB(a:) = (1 — x)9{l — x), where x = jk"^. For 

spinors this corresponds to a shape function rp such that x(l + rB(a:)) = x(l + rp(x))2. For it, 
the threshold functions can be computed analytically, and give 


I'rw) =. 


9 1 _!Z^ 

^ ^ d+2 


o 1 

,(F)d. ^ ^ ^ 


d 1 + 00 ’ 

J]±_ 
d+1 






d 


m. 




ni 


1 


1 _ r-i/i 

d+1 


+ n2- 


1 _ dl±. 

^ d+2 


(1 + ooiy+^^ (1 + 002)'^^ (1 + ooi)^^ (1 + a;2)^+^2 


m 


(F)d, 


to) = 


to) = 


1 


+ 


1 -??V’ 


{l+ooY {d-2){l + uof 

^ ^ (1+Cdl)4’ 

1 _ Iht. 

(FB)d/ N d+1 

”*'4 (“ 1 -“^) = (l+^.)(l+^,) 2 - 


1 - hi/- 1 

2d-4 ^ 47 (1 +w)2 ’ 
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Appendix B: Properties of the large-A^: solutions 


In both versions of the LPA, with or without = 0, and also in the LPA', Eq. (III. 8 ) enables 
us to write the potentials in the form 


4'yd 


h{4>) =Ch(l)^ , v{4>) = -] 


d . d 


„2 i2n 


(B.l) 


The behavior of v for (p —>■ ±oo is 


?;asympt(</>) ^ IPf'' (B.2) 

and, since we are assuming 2 < d < 4, the gamma function in front of |c/i|'^ is positive. If c„ ^ 0, 
the scalar potential can be real only if (—1)'^’^ has a real branch, that is if 

TTL 

d = — , j G {1,3, 5, ... } , m G N , 2nj <m < Anj . (B.3) 

nj 

Its stability further requires 

^d-\-l—d/2 

\ch\'^ - r(-d/2) (B.4) 

and for special values of Cy and Ch, namely when \ch\ = |c/i^crit|, it can become asymptotically 
flat (possibly only on one side) instead of growing like 0'^”. 

In order to understand the physical properties of the large-A/- FP’s, we need to consider 
the RG flow in vicinity of the corresponding critical points. In particular we consider the 
linearization of the flow, by looking at small fluctuations of the potentials v = v + 5v, h = h + 6h 
and for eigenvalue solutions 


5v = —95v , 6h = —9Sh . 


These equations at large-Aj are extremely simple and, for the linearized regulator, they read 

(B.5) 
(B. 6 ) 


- 96v = -dSv + -pdv' + ^ 

n 2 d (1 + h?)^ 

-96h = -6h + -pdh' + . 

n 2 


In this appendix we want to sketch a study of the properties of these FPs as well as of the 
linearized flow around them. We believe it can be instructive to consider separately the results 
obtained with or without the inclusion of the flow equation for rj^. This will make evident that 
larger truncations, out of the reach of the present work, are necessary to get a complete picture 
of the large-Aj multicritical Yukawa theories. 
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v(0) 



FIG. 14: The d = 3, n = 2, FP scalar potential at nonvanishing c„. Left panel: c„ = — 1 and Ch € 
{c?t,crit + 10“^, c?t_crit, c?t,crit — 10“^}, from bounded (green) to unbounded (blue). Right panel: = 1 

and Ch G {ch, crit + 2, Cft,,crit, c?i,crit - 2}, from steeper (green) to broader (blue). 



FIG. 15: The d = 8/3, n = 3, FP scalar potential at nonvanishing c„. Left panel: Cy = —1 and 
Ch G {ch,crit + 10~^, Ch,crit, Ch,crit — 10“^}, from bounded (green) to unbounded (blue). Right panel: 
Cy = 1 and Ch G {ch,crit + 2, Ch,crit, Ch,crit - 2}, from steeper (green) to broader (blue). 


1. LPA 


If we set by hand = 0, regardless of c„ or Ch Eq. (III.8) leaves a discrete set of dimensions 


as the only possibility, the ones in Eq. (III.ll). As a consequence dn = 2(n + 1) and the scalar 
potential is real and even also in case Cy / 0. The stability properties, depending on Ch and 


Cy according to Eq. (B.4), are illustrated in the plots of Fig. 14 and Fig. 15 The special case 


= 0 is shown in Fig. 16 


Let us now turn to the linear perturbations of these FP’s. By definition in the LPA one 
neglects a possible change of anomalous dimension. Thus, setting 6r]^ = 0, the solution for the 
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FIG. 16: The even FP scalar potentials for c„ = 0. For illustration the value of Ch has been chosen 


according to Eq. (B.12|, even if this is mandatory only for n = 1 in the LPA'. Left panel: n = 1 and 
d £ {3.5,3, 2.5}, from steeper (green) to broader (blue). Right panel: n £ {2,3,4, 5}, in the corresponding 
dimension d = 2 + 2/n, from steeper (green) to broader (blue). 


perturbations reads 


Sh((^) = 5ch(j) 


N 


(B.7) 


6v{cP) = 


1 


d 




d-2 


d 


d 


2 i2n 


fa 1 , 1---2 


. 1 

Here we restricted our analysis to the perturbations with 6ch A 0) required their smoothness 
by setting (1 — Q)n = N £ N. For the special case 6ch = 0 the solution is simply Sv{4>) = 6cv(j)^ 
with critical exponent 6m = d — M/n, and will not be discussed any further. Notice that these 
eigenfunctions are independent of c„, which is due to the suppression of scalar loops in the 
laige-Xf limit. They are regular at the origin, since the leading behavior is 


Recall that the FP potential had, as leading small field dependence, as a consequence, 
the relevant perturbations with N < n change the behavior of the potential at the origin, the 
marginal ones only change the coefficient in front of and the irrelevant ones leave the leading 
term unaltered. For large value of the field 


Sv(4.) (&„ + sgn(^)*> (B.9) 

while the FP potential behaves like at infinity. As a consequence, the irrelevant perturba¬ 
tions with N > n completely change the asymptotic behavior of the potential for large fields, 
the marginal ones with N = n only change the coefficient in front of the leading power, and the 
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FIG. 17: The d = 3 and n = 1 FP scalar potential at nonvanishing Cy. Left panel: Cy G 
{c„,crit ~ 1; c„_crit; c„_crit + 1}, from bounded (green) to unbounded (blue). Right panel: c„ G {—Cy^crit — 
Ij —c„,crit) —Cy^crit +1}, from bounded (green) to unbounded (blue). Notice that the value of the potential 
at the origin is arbitrary, while its behavior for large helds is not. 


relevant ones only change the sub-leading terms. Clearly this is not the case for those potentials, 
with special values of c„, that are asymptotically flat. 

Let’s now discuss the symmetry properties of the perturbations. Trivially, the symmetry of 
Yukawa potential under Z 2 is preserved or violated depending on n and N. We now want to 
understand what this entails for the scalar potential. Recall that in the LPA dn = 2(n -|- 1) 
and the FP v is always even. Then, the fluctuations behave as and whenever N + n 

is odd, the Z 2 symmetry of both h and v at the FP is spoiled by the perturbations. Among 
these symmetry breaking perturbations, the irrelevant ones, with N > n, give rise to unstable 
potentials. Notice that the relevant perturbations, even if spoiling symmetry, do not directly 
cause instabilities (though they might induce them indirectly, i.e. beyond linearization). The 
possibility to have stable theories with no definite Z 2 symmetry emanating from symmetric FPs 
in the UV or IR is in any case a question that requires a global study of the RG flow, and it is 
beyond the scope of this work. 


2. LPA' 


So far we have not used the flow equation for In order to do so, we first have to analyze 
the possible presence of a nontrivial minimum for v. The general expectation is that, since only 
fermion loops survive in the leading order of the 1/Xf expansion, the potential is always in the 
symmetric regime. This is suggested by the expansion of the potential around the origin, based 


on Eq. (III.7). We assume that this is always the case for the time being, as it is indeed for 
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FIG. 18: The d = 7/3 and n = 1 FP scalar potential at nonvanishing Cy. Left panel: c„ G 
{c„,crit — Ij c„_critj c„_crit + 1}, from bounded (green) to unbounded (blue). Right panel: Cy G {—Cy^crit — 
1) —Cy^crit! —Cy,crit + 1}) from bounded (green) to unbounded (blue). 


every specific example we have considered. Under this assumption, we need to take the (/>—)■ 0 
limit of the equation for which is proportional to /i'((())^, i.e. to Therefore, only 

for n = 1 such a limit can be nonvanishing. This shows how LPA' is an improvement of LPA 
only for the n = 1 critical theory. For the remaining values of n, one finds again rj^ = 0, which 
artificially forces the dimension d to its critical value. We expect this condition to be lifted by 
more general truncations, and a nontrivial should emerge for any n. 

Let’s then discuss the change brought by LPA' in the description of the large-Aj n = 1 FP. 


As argued in Sect. Ill the nontrivial allows for the existence of the non GauBian FP in any 
d < 4, as long as 


ry^ = 4- d, n = l. 


(B.IO) 


Actually this is the case only for the Z 2 symmetric solution with Cy = 0. As soon as 7 ^ 0 the 
reality of the potential requires 


m 


d = — , j E {1, 3, 5,... } , m E N , 2j < m < Aj . 
J 


(B.ll) 


Regardless of Cy, by using Eq. (B.IO) the flow equation for can be solved for Ch as a function 
of d, giving [16] 


2 _ d(4-d)(d-2) 
— 


Vd{6d - 8 ) 


(B.12) 


Then, the stability condition Eq. (B.4) for the nonvanishing Cy EP’s is best phrased as a bound 
on Cy 


Cy ^ ) Cy^crit — 


r(-d/ 2 ) 

2d+lj^d/2 


d(4-d)(d-2) 

Vd{6d - 8 ) 


d/2 


(B.13) 
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and additionally, only for odd m, < c^,crit- The scalar FP potential with / 0 is an even 
function if and only if m is even. 

Let us then turn to perturbations, and allow for a nontrivial br]^. We postpone for a while 
the task of solving the linearized equation for which provides us the first correction brj^ to 
the anomalous dimension, as a function of the FP h and 5h. This is because such an equation 
involves the variation 5(/>o in the location of the minimum of the potential, which in turn can be 
computed from the variation of the potential by the formula 


b^Q 


(5u'(0) 

u"(0) 


(B.14) 


where we stuck to our assumption that the minimum of the FP potential is always trivial. As a 


consequence we hrst solve for bv and bh as parametric functions of brj^, and then plug Eq. (B.14) 
into the linearized equation for r^j), to compute the actual br](j). Solving for bh is again trivial, 
and it immediately allows us to extract the eigenvalues of the linearized flow. When 9^0 the 
solution for bh is 


bh{(j)) = bch4> - 


% 


n 


2 n — N 


Ch(j)^ , iV G N , N^n 


(B.15) 


where again we focused on bch 7 ^ 0 and set iV = (1 — 9)n G N. For 0 = 0 instead 

bh{(l>) = bchcf)"' - log((/)) 


(B.16) 


Notice that the second term in the last equation is simply the first order in the expansion of 
which is the exactly marginal /i, around the n-th FP. As a consequence, the 


apparent instability that can come from the second term in Eq. (B.16) is actually a fake of 
linearization, as long as brj^ > —^jn. On the other hand, a logarithmic singularity at the 
origin appears even beyond linearization, and we believe this to be a pathology produced by 
the leading order in XjXj. The solution to this pathology will come soon, in the form of the 
constraint brjfj, = 0 for these perturbations. 

The equation for bv is much more involved in the LPA' than in the LPA, since it now depends 
on the EP potential. Yet, its solutions for generic bij^ can be easily given analytically. It is not 
necessary to show them here. It suffices to report that quite in general they have the property 
bv'{0) = 0 , as it could be expected by the argument that fermion loops are generally associated 
with scalar potentials with a trivial minimum As a consequence the scalar potential stays in 


For 5ch 7 ^ 0 only the n = 1, = 0 case gives rise to a nonvanishing h«'(0). For Sch = 0 only the M — 1 case. 

In what follows we discard these cases. 
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the symmetric regime. Notice that this does not entail that the 5v{4>) is also in the symmetric 
regime. 

With this piece of information, one can work out the linearized 6r]^, by varying the r.h.s. of 
the flow equation for with respect to h and v (whose fluctuations still depend parametri¬ 
cally on 6r]^ itself) and while keeping cpQ fixed, and then taking tpQ —)> 0. The latter limit 
makes the r.h.s. vanishing unless n = 1, in which case it reaches a d-dependent constant times 
c\5r]fj,. ^ Hence, for general n and N we find = 0, which boils the analysis of the linearized 
perturbations down to the one sketched in the last Section within the LPA. 


3. d = 4 


The expression in Eq. (B.l), cannot be used in d = 4 nor in d = 2, since the hypergeometric 


function in v has simple poles at these values. The d —>■ 2 case is out of the reach of the present 
paper. In the d —?• 4 limit, instead, the canonical dimensional terms survive also in the LPA, and 
by integrating the large-Aj system of flow equations one can find the following FP solutions 

h{(p) = Ch4>'^ , v{(t)) = Cv log(c| -h (B.17) 


where we already demanded the Yukawa potential to be smooth, according to Eq. (III.8). The 
crucial fact is again that the minimum of v is always trivial. This allow us to take the lyio 0 
limit of the equation for r]^. For n = 1 this leaves us with the equation (^ = il(j> = 0) thus 
boiling every feature of the critical theory down to the classical counting. For n > 2 we hnd the 


constraint = 0, which is inconsistent with Eq. (III.8) and therefore eliminates these solutions 


[1] D. J. Gross and A. Neveu, Phys. Rev. D 10 (1974) 3235. 

[2] B. Rosenstein, B. J. Warr and S. H. Park, Phys. Rev. D 39 (1989) 3088; A. Kovner, B. Rosenstein 
and G. Gat, Helv. Phys. Acta 65 (1992) 411. 

[3] B. Rosenstein, H. L. Yu and A. Kovner, Phys. Lett. B 314 (1993) 381. 

[4] J. Zinn-Justin, Nucl. Phys. B 367 (1991) 105; J. Zinn-Justin, “Phase transitions and renormalization 
group”, Oxford, UK: Oxford Univ. Pr. (2007) 452 p 

® Such a constant is actually infinite for the marginal perturbation, the r.h.s. inheriting a logarithmic singularity 
at the origin from Sh. Yet the simple way to cure this pathology and get a self-consistent answer is to set 
Srj^ = 0 . 


49 








[5] J. A. Gracey, Nucl. Phys. B 341 (1990) 403; Nucl. Phys. B 367 (1991) 657; C. Luperini and P. Rossi, 
Annals Phys. 212 (1991) 371. 

[6] A. N. Vasiliev and M. I. Vyazovsky, Theor. Math. Phys. 113 (1997) 1277 [Teor. Mat. Fiz. 113 (1997) 
85]; J. A. Gracey, Nucl. Phys. B 802 (2008) 330 [arXiv:0804.1241 [hep-th]]. 

[7] J. A. Gracey, Int. J. Mod. Phys. A 9 (1994) 567 [hep-th/9306106]. 

[8] J. A. Gracey, Int. J. Mod. Phys. A 9 (1994) 727 [hep-th/9306107]; A. N. Vasiliev, S. E. Derkachov, 

N. A. Kivel and A. S. Stepanenko, Theor. Math. Phys. 94 (1993) 127 [Teor. Mat. Fiz. 94 (1993) 
179]. 

[9] S. Hands, A. Kocic and J. B. Kogut, Annals Phys. 224 (1993) 29 [hep-lat/9208022]. 

[10] L. Karkkainen, R. Lacaze, P. Lacock and B. Petersson, Nucl. Phys. B 415 (1994) 781 [Erratum-ibid. 
B 438 (1995) 650] [hep-lat/9310020]. 

[11] E. Focht, J. Jersak and J. Paul, Nucl. Phys. Proc. Suppl. 47 (1996) 709 [hep-lat/9509040]; S. ChristoH 
and G. Strouthos, JHEP 0705 (2007) 088 [hep-lat/0612031]. 

[12] S. Ghandrasekharan and A. Li, Phys. Rev. D 88 (2013) 021701 [arXiv:1304.7761 [hep-lat]]. 

[13] L. Wang, P. Gorboz and M. Troyer, New J. Phys. 16 (2014) 10, 103008 [arXiv: 1407.0029 [cond- 

mat.str-el]]. 

[14] L. Rosa, P. Vitale and G. Wetterich, Phys. Rev. Lett. 86 (2001) 958 [hep-th/0007093]. 

[15] F. Hofling, C. Nowak and C. Wetterich, Phys. Rev. B 66 (2002) 205111 [cond-mat/0203588]. 

[16] J. Braun, H. Gies and D. D. Scherer, Phys. Rev. D 83 (2011) 085012 [arXiv:1011.1456 [hep-th]]; J. 
Phys. A 46 (2013) 285002 [arXiv: 1212.4624 [hep-ph]]. 

[17] H. Sonoda, Prog. Theor. Phys. 126 (2011) 57 [arXiv:1102.3974 [hep-th]]. 

[18] L. Janssen and 1. F. Herbut, Phys. Rev. B 89 (2014) 20, 205403 [arXiv: 1402.6277 [cond-mat.str-el]]. 

[19] J. Borchardt and B. Knorr, arXiv:1502.07511 [hep-th]. 

[20] N. Dorey and N. E. Mavromatos, Phys. Lett. B 250 (1990) 107; 1. J. R. Aitchison and N. E. Mavro- 
matos, Phys. Rev. B 53 (1996) 9321 [hep-th/9510058]; N. E. Mavromatos and J. Papavassiliou, Re¬ 
cent Res. Devel. Phys. 5 (2004) 369 [cond-mat/0311421]. M. Franz and Z. Tesanovic, Phys. Rev. Lett. 
87 (2001) 257003; 1. F. Herbut, Phys. Rev. Lett. 88 (2002) 047006 [cond-mat/0110188]; Phys. Rev. 
B 66 (2002) 094504 [cond-mat/0202491]; Phys. Rev. Lett. 94 (2005) 237001 [cond-mat/0410557]. 

[21] G. W. Semenoff, Phys. Rev. Lett. 53 (1984) 2449; 1. F. Herbut, Phys. Rev. Lett. 97 (2006) 
146401 [cond-mat/0606195]; 1. F. Herbut, V. Juricic and B. Roy, Phys. Rev. B 79 (2009) 085116 
[arXiv:0811.0610 [cond-mat.str-el]]; V. P. Gusynin, S. G. Sharapov and J. P. Garbotte, Int. J. Mod. 
Phys. B 21 (2007) 4611 [arXiv:0706.3016 [cond-mat.mes-hall]]; S. Hands and C. Strouthos, Phys. 
Rev. B 78 (2008) 165423 [arXiv:0806.4877 [cond-mat.str-el]]; J. E. Drut and T. A. Lahde, Phys. 
Rev. Lett. 102 (2009) 026802 [arXiv:0807.0834 [cond-mat.str-el]]; Y. Araki and T. Hatsuda, Phys. 
Rev. B 82 (2010) 121403 [arXiv: 1003.1769 [cond-mat.str-el]]. 

[22] H. Gies and L. Janssen, Phys. Rev. D 82 (2010) 085018 [arXiv:1006.3747 [hep-th]]; L. Janssen and 


50 



H. Gies, Phys. Rev. D 86 (2012) 105007 [arXiv:1208.3327 [hep-th]]. 

[23] T. Grover, D. N. Sheng and A. Vishwanath, Science 344 (2014) 6181, 280 [arXiv: 1301.7449 [cond- 
mat.str-el]]. 

[24] L. Balents, M. P. A. Fisher and G. Nayak, Int. J. Mod. Phys. B 12 (1998) 1033 [cond-mat/9803086 
[cond-mat.supr-con]], S. S. Lee, Phys. Rev. B 76 (2007) 075103 [cond-mat/0611658]. 

[25] Y. Nambu and G. Jona-Lasinio, Phys. Rev. 122 (1961) 345; J. B. Kogut, M. A. Stephanov and 
C. G. Strouthos, Phys. Rev. D 58 (1998) 096001 [hep-lat/9805023]. S. Chandrasekharan, J. Cox, 
K. Holland and U. J. Wiese, Nucl. Phys. B 576 (2000) 481 [hep-lat/9906021]; S. Hands, B. Lucini 
and S. Morrison, Phys. Rev. D 65 (2002) 036004 [hep-lat/0109001]. 

[26] D. U. Jungnickel and C. Wetterich, Phys. Rev. D 53 (1996) 5142 [hep-ph/9505267]. 

[27] H. Gies and M. M. Scherer, Eur. Phys. J. C 66 (2010) 387 [arXiv:0901.2459 [hep-th]]. 

[28] H. Gies, S. Rechenberger and M. M. Scherer, Eur. Phys. J. C 66 (2010) 403 [arXiv:0907.0327 [hep- 
th]]; H. Gies, L. Janssen, S. Rechenberger and M. M. Scherer, Phys. Rev. D 81 (2010) 025009 
[arXiv:0910.0764 [hep-th]]; H. Gies, S. Rechenberger, M. M. Scherer and L. Zambelli, Eur. Phys. J. 
C 73 (2013) 2652 [arXiv: 1306.6508 [hep-th]]. 

[29] Y. Nambu, in: Z. Ajduk, S. Pokorski, A. Trautman ed. (World Scientific Pub., ISBN 9971-50-691- 
2, 1988, p. 1-10) 11. Warsaw symposium on elementary particle physics: new theories in physics, 
Kazimierz (Poland), 1988; V. A. Miransky, M. Tanabashi and K. Yamawaki, Phys. Lett. B 221 
(1989) 177; Mod. Phys. Lett. A 4 (1989) 1043; H. Gies, J. Jaeckel and C. Wetterich, Phys. Rev. D 
69 (2004) 105008 [hep-ph/0312034]. 

[30] O. Zanusso, L. Zambelli, G. P. Vacca and R. Percacci, Phys. Lett. B 689 (2010) 90 [arXiv:0904.0938 
[hep-th]]; G. P. Vacca and O. Zanusso, Phys. Rev. Lett. 105, 231601 (2010) [arXiv: 1009.1735 [hep- 
th]]. 

[31] J. M. Pawlowski and F. Rennecke, arXiv:1403.1179 [hep-ph]; J. Braun, L. Fister, J. M. Pawlowski 
and F. Rennecke, arXiv:1412.1045 [hep-ph]. 

[32] C. Wetterich, Phys. Lett. B 301, 90 (1993). 

[33] K. Aoki, Int. J. Mod. Phys. B 14, 1249 (2000); J. Braun, J. Phys. G 39, 033001 (2012) 
[arXiv:1108.4449]; J. Berges, N. Tetradis and C. Wetterich, Phys. Rept. 363, 223 (2002) [hep- 
ph/0005122]; B. Delamotte, Lect. Notes Phys. 852, 49 (2012) [cond-mat/0702365 [COND-MAT]]; 
P. Kopietz, L. Bartosch and F. Schutz, Lect. Notes Phys. 798, 1 (2010); T. R. Morris, Prog. 
Theor. Phys. Suppl. 131 (1998) 395 [hep-th/9802039]; O. J. Rosten, Phys. Rept. 511 (2012) 177 
[arXiv:1003.1366 [hep-th]]. 

[34] G. Felder, Commun. Math. Phys. Ill, 101-121 (1987). 

[35] A. Codello, J. Phys. A 45 (2012) 465006 [arXiv: 1204.3877 [hep-th]]; A. Codello and G. D’Odorico, 
Phys. Rev. Lett. 110 (2013) 141601 [arXiv: 1210.4037 [hep-th]]; A. Codello, N. Defenu and 
G. D’Odorico, arXiv: 1410.3308 [hep-th]. 


51 



[36] T. R. Morris, Phys. Lett. B 334 (1994) 355 [hep-th/9405190]. 

[37] T. R. Morris, Phys. Lett. B 345 (1995) 139 [hep-th/9410141]. 

[38] T. R. Morris, Phys. Lett. B 329 (1994) 241 [hep-ph/9403340]; T. R. Morris and M. D. Turner, Nucl. 
Phys. B 509 (1998) 637 [hep-th/9704202]. 

[39] F. Synatschke, J. Braun and A. Wipf, Phys. Rev. D 81 (2010) 125001 [arXiv: 1001.2399 [hep-th]]. 

[40] H. Gies, F. Synatschke and A. Wipf, Phys. Rev. D 80 (2009) 101701 [arXiv:0906.5492 [hep-th]]. 

[41] M. Heilmann, T. Hellwig, B. Knorr, M. Ansorg and A. Wipf, JHEP 1502 (2015) 109 [arXiv: 1409.5650 
[hep-th]]. 

[42] T. Hellwig, A. Wipf, and O. Zanusso, in preparation. 

[43] D. Bashkirov, arXiv:1310.8255 [hep-th]. 

[44] O. Aharony, A. Hanany, K. A. Intriligator, N. Seiberg and M. J. Strassler, Nucl. Phys. B 499 (1997) 
67 [hep-th/9703110]. 

[45] D. Mesterhazy, J. Berges and L. von Smekal, Phys. Rev. B 86 (2012) 245431 [arXiv: 1207.4054 
[cond-mat.str-el]]. 

[46] R. Percacci and G. P. Vacca, arXiv: 1501.00888 [hep-th]. 


52 



